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ABSTRACT 


There  is  a  great  need  for  a  large  and  accurate  global  data  base  in  the  critical  60  to  120  km  region  of 
the  atmosphere  in  order  to  understand  in  detail  the  spatial  and  temporal  variabilities  in  its  mass 
density.  The  only  practical  method  of  providing  global  information  is  to  use  a  remote  sensing 
technique  from  earth  orbit.  Several  possible  techniques  had  been  proposed  which  have  the 
potential  to  measure  the  density  in  this  important  atmospheric  region.  These  include  measuring  the 
brightness  of  sunlight  that  is  rayleigh  scattered  by  the  molecular  nitrogen  and  oxygen,  measuring 
the  absorption  of  starlight  during  stellar  occultation,  and  measuring  the  intensity  of  molecular 
oxygen  resonance  fluorescence.  A  feasibility  study,  funded  by  the  Air  Force  Geophysics 
Laboratory,  has  been  conducted  here  at  the  Space  Physics  Research  Laboratory,  University  of 
Michigan,  on  the  relative  merits  of  these  techniques.  This  report  presents  the  results  of  the  trade¬ 
off  studies  wherein  various  sources  of  error  and  complexity  of  the  data  analysis  are  evalulated.  An 
instrumental  conceptual  design  is  also  prepared  identifying  the  major  hardware  elements  required  to 
implement  the  remote  sensing  technique. 


1.  INTRODUCTION 


Orbital  reentry  of  large  maneuverable  vehicles  and  the  potential  to  modify  orbital  parameters  using 
the  aerobreaking  technique  have  increased  the  need  to  understand  in  detail  the  variability  of 
atmospheric  density  in  the  mesosphere  and  lower  thermosphere.  Acceleration  data  (Champion, 
1986)  obtained  during  shuttle  reentry  illustrates  the  fact  that  the  atmosphere  exhibits  variation 
which  often  exceeds  expectation  and  which  raise  questions  regarding  the  safe  return  of  heavy 
Shuttle  payloads  from  polar  orbit. 

The  variability  of  mass  density  in  the  mesosphere  has  been  apparent  for  many  years,  being  seen  in 
sounding  rocket  data,  backscatter  radar  profiles,  and  airglow  optical  emissions  (Wasser  and 
Donahue,  1979).  The  systematic  variations  that  have  been  determined  from  these  data  sets  are  the 
obvious  seasonal,  latitudinal,  and  local  time  dependencies  that  are  prescribed  in  the  major  models 
in  use  today.  These  models,  however,  are  based  on  limited  data  and  thus,  are  not  useful  in  making 
predictions  or  forecasts  of  the  effects  of  planetary  waves  or  gravity  wave  breaking.  The  wave 
processes  will  never  be  accurately  predicted,  but  like  the  weather  in  the  troposphere  we  may 
ultimately  be  able  to  base  a  forecast  on  operational  information.  The  primary  factor  holding  back 
the  improvement  of  these  models  and  the  development  of  a  forecasting  technique  is  the  lack  of  a 
sufficiently  large  and  accurate  morphological  data  base. 

This  global  data  base  must  rely  on  a  sampling  technique  that  will  provide  a  morphology  of  the 
atmospheric  density  in  the  critical  60  to  120  km  region,  with  the  consistency  required  to  quantify 
the  higher  frequency  variability.  The  only  practical  method  of  providing  global  information  is  to 
use  a  remote  sensing  technique  from  earth  orbit,  given  a  technique  that  has  the  required  accuracy 
and  spacial  sensitivity.  Such  a  technique  could  be  tested  using  first  sub-orbital  evaluation, 
followed  by  shuttle  validation,  and  finally  deployed  on  a  satellite  in  earth  orbit  to  provide  the 
required  global  information. 

Several  possible  remote  sensing  techniques  have  been  proposed  which  have  the  potential  to 
measure  the  atmospheric  density  in  this  important  atmospheric  region.  These  include  measuring 
the  brightness  of  sunlight  that  is  rayleigh  scattered  by  the  molecular  nitrogen  and  oxygen  (Rusch  et 
al„  1983),  measuring  the  absorption  of  starlight  during  stellar  occultation  (Hays  and  Roble,  1968). 
and  a  new  technique  where  the  fluorescence  of  molecular  oxygen  is  measured  photometrically. 
Each  of  these  methods  of  sensing  the  density  has  inherent  advantages  as  well  as  weaknesses.  In 
this  report,  we  will  examine  these  techniques  in  details  by  simulating  the  measurements. 
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developing  the  density  recovering  algorithms,  designing  conceptually  the  best  suitable  instrument, 
and  analyzing  the  errors  in  the  recovered  density.  The  three  proposed  approaches  will  be 
evaluated,  focusing  primarily  on  the  accuracy  of  the  measurements,  the  complexity  of  the  data 
processing  and  analysis. 
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2  .  STELLAR  OCCULT ATION  TECHNIQUE 


Stellar  occultation  has  been  used  for  many  years  to  determine  the  composition  and  number  density 
profiles  in  the  Earth  and  planetary  atmospheres.  During  the  early  years,  most  of  the  occultation 
measurements  have  been  made  from  rockets  by  examining  the  absorption  of  solar  ultraviolet 
radiation  in  certain  atmospheric  bands.  For  example,  vertical  ozone  density  profiles  were  obtained 
by  Johnson  et  al.  (1951)  from  series  of  rockets  using  this  technique.  Later,  artificial  satellites  were 
used  (Venkatestwaran  et  al.  1961;  Hinteregger  and  Hall,  1969;  Hays  and  Roble,  1973). 

Basically,  the  occultation  technique  is  similar  to  the  classical  technique  of  absorption  spectroscopy, 
and  has  been  investigated  in  detail,  with  the  early  theoretical  studies  of  Hays  and  Roble  (1968a), 
examination  of  errors  by  Roble  and  Hays  (1972),  and  recent  reviews  by  Atreya  (1981).  As  shown 
in  Fig.  2.1,  the  Sun,  Moon,  or  a  star  is  used  as  a  source  of  light  and  the  satellite  photometer  as  a 
detector,  with  the  atmosphere  between  acting  as  the  extinction  cell.  At  long  wavelengths  (such  as 
visible  and  radio  wave)  the  extinction  results  primarily  from  differential  refraction  of  the  light  by 
atmospheric  particles.  At  shorter  wavelengths  (e.g.  ultraviolet),  however,  absorption  by  the 
various  atmospheric  species  is  mainly  responsible  for  the  extinction  of  light.  As  the  satellite  moves 
in  its  orbit,  the  intervening  atmosphere  extincts  progressively  more  of  the  light  as  it  traverses 
deeper  layers  of  the  atmosphere.  The  occultation  intensity  data  can  be  related  to  the  tangential 
column  number  density  of  the  atmospheric  spieces  provided  that  the  extinction  cross-sections  are 
known.  The  equation  for  the  tangential  column  number  density  can  then  be  easily  inverted  to  give 
the  number  density  profile  at  the  occultation  tangent  point.  It  is  worthwhile  to  note  that  the 
atmospheric  densities  obtained  from  this  technique  are  independent  of  the  knowledge  of  absolute 
flux  of  the  light  source,  as  the  technique  depends  only  on  the  relative  attenuation  of  the  light. 
However,  the  nature  of  the  light  source  does  affect  the  design  of  the  instrument  in  terms  of  the 
accuracy  required  for  the  recovered  density  profile. 

In  this  chapter  we  describe  a  technique  for  retrieving  the  temperature  and  the  molecular  oxygen 
number  density  profiles  in  the  critical  60  to  120  km  altitude  region  from  occultation  intensity  data. 
A  data-reduction  technique  similar  to  the  one  described  by  Hays  and  Roble  (1968a)  is  developed, 
and  a  theoretical  error  analysis  is  used  to  determine  the  accuracy  of  the  retrieved  temperature  and 
density  profiles  in  the  presence  of  random  statistical  error  in  the  measured  intensities.  This  allows 
us  to  examine  the  retrieval  accuracy  as  a  function  of  the  star  brightness  and  its  color  type. 
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2.1  Principle  of  Stellar  Occultation 

A  detector  on  board  a  spacecraft,  as  shown  in  Fig.  2.1,  measures  the  unattenuated  intensity,  loo, 
of  the  source  when  the  absorption  path,  the  line  connecting  the  source  and  the  detector,  is  above 
the  atmosphere.  As  the  spacecraft  moves  in  its  orbit,  the  absorption  path  descends  deeper  into  the 
atmosphere,  resulting  in  progressively  greater  absorption  of  light.  The  attenuation  of  mono¬ 
chromatic  light  of  wavenumber  v  due  to  absorption  by  the  i1*1  atmospheric  species  is  given  by 
Beer's  Law,  which  states 

I  (v,  rc)  =  IM  (v)  exp  [  -  x.  (v)  ],  (2.1) 


where  I  (v,r0)  is  the  attenuated  intensity  at  tangent  point  of  radius  r0,  and  x*  (v)  is  the  line  of  sight 
optical  depth  of  species,  i.  If  we  assume  that  the  absorption  cross-section  of  the  i*  atmospheric 
species  at  wavenumber  v,  O;  (v),  is  independent  of  altitude,  x,  (v)  can  be  written  as 

x.  (v)  =  a  (v)  N.  fro),  (2.2) 

where  N;  (r0)  is  its  line  of  sight  column  density  of  the  absorber  N;  at  radius  rD.  In  the  case  of  a 
spherically  stratified  atmosphere,  the  tangential  column  density  N,  (rQ)  is 


n.  (r)  r  dr 


(2.3) 


Figure  2.1 
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Figure  2 . 2 


where  n,(r)  is  the  number  density  of  the  i!h  absorbing  species  at  a  distance  r  from  the  center  of  the 
Earth.  This  equation  is  the  Abel  integral  equation  (c.f.  Hays  and  Roble,  1968  a,  b)  which  can  then 
be  directly  inverted  to  give  the  number  density  of  the  absorbing  species  at  radius  r 


n,r.  If  idN(r0)/dr0  ]  dr0 


Hays  and  Roble  (1968b)  calculated  the  tangential  ultraviolet  transmission  for  the  mesosphere  and 
the  lower  thermosphere.  They  showed  that  in  certain  specific  wavelength  and  altitude  intervals  the 
ultraviolet  starlight  is  absorbed  by  a  single  species  and  attenuations  due  to  refraction  and  Rayleigh 
scattering  are  negligible.  These  include  a  region  near  1500A  (66000  cm1)  where  molecular 
oxygen  acts  as  a  single  absorber  in  the  lower  thermosphere  and  also  a  region  near  25CK)A  (40000 
cm1)  where  ozone  is  a  single  absorber  in  the  mesosphere. 


Among  various  photoabsorption  processes  by  the  molecular  oxygen,  photodissociations  through 
the  Herzberg  continuum,  Schumann-Runge  bands  and  Schumann-Runge  continuum  make  the 
most  contributions  to  the  extinction  of  solar  ultraviolet  light.  Fig.  2.2  shows  the  required 
magnitudes  of  absorption  cross-section  as  a  function  of  tangent  height  which  would  give  optical 
depths  of  0.5,  1.0  and  1.5  for  an  O2  atmosphere  during  an  occultation  sequence.  The  absorption 
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cross-sections  of  molecular  oxygen  due  to  the  three  processes  are  also  shown.  The  observations 
near  the  Schumann-Runge  continuum  allow  measurement  of  density  down  to  about  120  km.  This 
has  been  successfully  demonstrated  in  the  observations  of  the  OAO-A2  satellite  reported  by  Hays 
and  Roble  (1973).  As  the  altitude  decreases  and  the  tangential  O2  column  density  increases,  the 
spectral  regions  with  smaller  absorption  cross-sections  are  required,  and  the  various  Schumann- 
Runge  bands  between  1750A  and  2000A  (50000-57000  cm1)  and  the  Herzberg  continuum  beyond 
2000A  (<  57000  cm'1)  become  applicable.  Although  these  regions  were  also  observed  by  the 
OAO-A2  satellite,  the  occultation  data  were  not  analyzed  due  to  the  complexity  of  the  band  system 
and  the  uncertainty  of  knowledge  of  its  relatively  small  cross-sections.  Recently  cross-sections  of 
the  Herzberg  continuum  (Cheung  et  al.,  1986a)  and  the  Schumann-Runge  bands  as  well  as  the 
positions  of  each  of  the  individual  rotational  lines  have  been  obtained  at  very  high  resolution 
(-0.04  cnv1)  from  photoabsorption  measurements  using  a  photoelectric  scanning  spectrometer 
(Yoshino  et  al.,  1983,  1984;  Cheung  et  a!.,  1984;  Yoshino  et  al.,  1987).  This  allows  us  to 
understand  the  band  systems  better  theoretically,  and  enables  us  to  examine  the  feasibility  of  using 
the  stellar  occultation  technique  to  extract  O2  number  density  in  the  critical  60  to  120  km  region. 


2.1.1  02  Absorption  Cross-sections 
Herzberg  Continuum 

The  Herzberg  continuum  is  one  of  the  two  important  dissociation  continua  in  02  photoabsorption 
processes.  It  occurs  near  the  O2  dissociation  limit  where  \<2420A  and  involves  a  forbidden 
transition  between  02(X3!  ,  ground  state  and  Q2(A3Zll  +  )  excited  state,  leading  to  photo¬ 
dissociation  products  of  two  0(3P)  atoms. 

The  absorption  at  this  transition  is  rather  weak,  thus  measurements  of  the  cross-sections  are  very 
difficult  (Buisson  et  al.,  1933;  Gotz  and  Maier-Leibnitz.  1933).  Ditchbum  and  Young  (1962)  used 
long  absorption  paths  to  increase  the  optical  depth  and  obtained  cross-sections  in  this  spectral 
region  varying  from  1.4xl0'24  cm2  at  2500A  to  2.64x10  23  cm2  at  2000A.  Recently,  using  a 
6.65m  photoelectric  scanning  spectrometer,  Cheung  et  al.,  (1986a)  were  able  to  measure  the 
Herzberg  continuum  cross-sections  more  accurately  in  the  wavelength  region  between  1950A  and 
2400A.  They  obtained  a  peak  cross-section  of  6.62X10'24  cm2  at  2010A,  which  is  smaller  than  the 
ones  that  Ditchbum  and  Young  (1962)  had  derived  and  Banks  and  Kockart  (1973)  had  adopted. 
In  this  study,  we  use  the  laboratory  measurement  of  Cheung  et  al.  (1986a). 


Schumann-Runge  Continuum 


The  stronger  Schumann-Runge  Continuum  arises  because  the  permitted  electric-dipole  transition 
02(X3Sg  )  +  hv  - >  O2  (B3Lu‘), 

for  wavelengths  shorter  than  1750A  down  to  1250A,  and  leads  to  a  final  molecular  state  lying 
above  the  threshold  for  the  dissociation 

02  (B3£u  ) - >  0(3p)  +  O('D). 

It  has  an  absorption  cross-section  that  peaks  with  a  value  of  1.1  x  1017  cm2  near  1400A  and 
decreases  toward  the  longer  wavelengths.  Early  measurements  by  Ladenburg  and  Van  Voorhis 
(1933)  and  by  Watanabe  et  al.  (1953)  differ  by  about  25  percent  at  some  wavelengths.  Later 
measurements  by  Metzger  and  Cook  (1964)  and  Huffman  et  al.  (1964)  also  differ  by 
approximately  10  percent.  Hall  et  al.  (1965)  then  indicated  the  need  for  more  accurate 
measurements,  given  the  differences  in  the  published  values  of  the  absorption  cross-sections,  and 
the  necessity  to  accurately  know  this  cross-section  for  stellar  occultation  measurements  of 
atmospheric  density.  In  an  attempt  to  investigate  the  effect  of  temperature  on  the  02  absorption 
coefficients,  Hudson  et  al.  (1966)  measured  the  absorption  cross-sections  of  the  Schumann-Runge 
Continuum  between  1580A  and  1950A  at  three  different  temperatures.  They  concluded  that  the 
maximum  change  in  the  absorption  cross-section  is  about  10%  when  the  temperature  of  the  oxygen 
is  increased  from  300  K  to  1000K.  In  our  study,  we  adopt  the  cross-sections  of  Metzger  and 
Cook  (1964)  for  1325A  <  X  <  1650A,  and  the  results  of  Hudson  et  al.  (1966)  at  300K  for  1650A 

<X<  1800A. 

Schumann-Runge  Bands 

The  Schumann-Runge  band  system,  originating  from  the  same  electronic  transition  as  the 
Schumann-Runge  Continuum,  begins  with  the  (0-0)  band  at  2026A  and  converges  to  the  second 
dissociation  limit  at  about  1750A.  Kreusler  (1901)  first  attempted  to  measure  the  absorption 
coefficient  in  this  spectral  region  and  obtained  a  cross-section  on  the  order  of  10'21  cm2.  Since  then 
these  bands  have  been  studied  extensively  by  several  investigators,  for  example,  by  Watanable  et 
al.  (1953)  and  by  Wilkinson  and  Mulliden  (1957).  Recently,  using  a  high  resolution  photoelectric 
scanning  spectrometer,  Yoshino  et  al.  (1983),  Cheung  et  al.  (1984)  and  Yoshino  et  al.  (1987)  were 
able  to  measure  the  absorption  cross-sections  of  each  of  the  individual  rotational  line  throughout 
the  wavelength  region  containing  the  (d',0)  bands  with  o'  =  0,12.  Errors  in  the  their  cross- 
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section  measurements  were  estimated  to  be  4  percent  and  somewhat  greater  for  cross-sections  less 
than  5.0  x  1022  cm2.  The  band  oscillator  strengths  were  obtained  at  two  temperatures,  300  K  and 
79K,  and  differences  between  them  were  observed  although  they  lie  within  the  measurement 
errors. 

Because  of  the  existence  of  a  repulsive  3nu  state,  whose  potential  curve  crosses  that  of  the  B3Zu  in 
the  vicinity  of  =  4,  the  lines  in  the  Schumann-Runge  bands  are  predissociative  lines,  i.e.  they  are 
unstable  against  molecular  dissociation.  As  a  result,  these  lines  are  broadened,  and  extensive 
overlaps  exist  between  rotational  lines  in  the  same  band  or  different  bands.  Determining  O2  density 
using  the  stellar  occultation  technique  from  the  photoabsorption  of  light  in  the  Schumann-Runge 
bands  thus  requires  line  by  line  calculations  based  on  a  band  model  that  will  take  into  account  the 
effect  of  temperature  on  both  the  band  oscillator  strength  and  the  absorption  coefficient  of  each 
rotational  line.  Next  we  will  describe  how  we  have  calculated  the  absorption  cross-section  of  each 
line. 

If  we  design;  -e  a  line  within  a  Schumann-Runge  band  by  its  upper  and  lower  states  of  the 
vibrational,  rotational,  and  total  angular  momentum  quantum  numbers,  o',  N',  J'  and  u",  N",  J”, 
respectively,  the  absorption  cross-section  O  (v)'\>'',N'N",J'J")  can  be  written  as 

a  (v'v",  N’N”,  IT)  =  ^  Rv2v.  S(V  v",  N’  N",  J'  J"),  (2.5) 

where  v  is  the  transition  frequency,  c  is  the  speed  of  the  light,  Ruv  is  the  transition  moment  which 
is  the  product  of  the  electronic  dipole  moment  and  the  vibrational  overlap  integral,  and 
S(i)V,N’N",JT’)  is  the  line  strength,  or  Hbnl-London  factor. 

The  oscillator  strengths  are  calculated  for  the  lines  in  the  Schumann-Runge  bands  arising  from  the 
F2  fine-structure  component  of  the  ground  X3Ig  state,  since  the  dependence  of  the  oscillator 
strength  on  the  fine-structure  component  is  expected  to  be  weak  (Lewis  et  al.,  1987).  The  Honl- 
London  factors  are  taken  from  Tatum  and  Watson  ( 1 97 1 )  for  a  coupling  case  intermediate  between 
Hund's  cases  (a)  and  (b).  They  are  calculated  using  the  accurate  O2  molecular  parameters  of  the 
X3£g-  state  (Steinbach  and  Gordy,  1975)  and  the  B3!^-  state  (Cheung  et  al.,  1986b).  All  of  the 
twelve  allowed  lines  in  the  P  and  R  branches,  including  the  weak  satellite  lines,  and  the  two 
forbidden  lines  corresponding  to  AN=±3  are  included.  The  vibrational  wave  functions  and  the 
overlap  integral  are  computed  based  on  the  X3Sg-  potential  of  Allison  et  al.  (1971)  and  the  B3Xu' 
potential  of  Allison  et  al.  (1986).  The  electronic  dipole  moment  of  Allison  et  al.  (1986)  was  also 
adopted. 
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The  shape  of  each  rotational  line  in  the  Schumann-Runge  bands  is  broadened  due  to 
predissociation  and  has  been  approximated  by  a  Lorentzian  profile  by  Hudson  and  Mahle  (1971) 
and  Omidvar  and  Frederick  (1987)  in  their  calculations  of  photodissiociation  rate.  The  absorption 
cross-section  of  a  line  i  at  wavenumber  v,  Oi(v),  using  the  Lorentz  distribution  function  is  given  by 


a  (v)  =  <r  (Vj) 


(AL/2)2 

(v  -  v.)2  +  (ALj/2)2  * 


(2.6) 


where  Vi  and  a(v,)  are  the  wavenumber  and  the  cross-section  at  the  line  center,  and  ALj  is  the 
FWHM  (Full  Width  at  Half  Maximum)  line  width.  The  validity  of  this  assumption  was  discussed 
in  detail  by  Omidvar  and  Frederick  (1987),  except  that  we  used  measured  Lorentz  widths  of 
Yoshino  (1988,  private  communication),  which  give  the  variation  of  AL,  with  respect  to  the  total 
angular  momentum  quantum  number  J. 

The  temperature  effect  on  the  absorption  coefficients  is  introduced  when  the  rotational  population 
of  02(X3rg-,\)"),  fN(T),  is  considered, 

fN  (T)  =  (2N  +  1)  exp  { -  E"/kT  }  /  Qr  (T),  (2.7) 

where  Qr  is  the  rotational  partition  function  given  by 

Qr  CO  *  T  (2N  +  1 )  exp  { -  E"/kT  } ,  (2.8) 

N=1A5 


T  is  the  rotational  temperature,  E"  is  the  energy  level  for  the  rotational  level  N,  and  N  takes  only 
odd  integer  numbers,  in  accordance  with  the  exclusion  principle  for  an  ordinary  molecule 

(Herzberg,  1950).  For  temperatures  below  300K  (the  case  in  the  lower  thermosphere  and 
mesosphere),  all  O2  molecules  are  populated  in  their  lowest  vibrational  level,  o"  =  0.  We  can  then 
define  an  effective  absorption  cross-section  ae,(v,  T)  as 

CTei  (v,  T)  =  fN  (T)  c.  (v)  (2.9) 

and  calculate  the  optical  depth  t(v)  by 

00  ao 

T(V)  =  II  n(s)  fN  (T)  a  (v)  ds  =  II  n(s)  Oei  (v,  T)  ds  (2.10) 

-00  *  1  -00 
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Figure  2.3 
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Figure  2.4 


to  consider  the  overlap  of  adjacent  lines.  All  of  the  cross-sections  appearing  hereafter  will  be 
referred  to  as  the  effective  cross-sections,  so  that  the  rotational  and  vibrational  populations  of  the 
C>2  molecules  do  not  have  to  be  considered. 

Fig.  2.3  shows  our  calculated  absorption  cross-sections  of  Schumann-Runge  (6,0)  band  in 
comparison  with  the  measurements  of  Yoshino  et  al.  (1987)  at  79K.  The  two  are  in  excellent 
agreement.  Fig.  2.4  presents  the  calculated  absorption  cross-sections  at  300K  for  the  same  (6,0) 
band.  As  the  temperature  increases,  the  lines  in  the  shorter  wavenumber  region,  representing  the 
transitions  from  higher  rotational  levels  of  the  ground  state  of  O2,  become  more  pronounced.  This 
apparent  change  in  band  structure,  which  appears  more  clearly  in  a  high  resolution  O2  absorption 
spectrum,  thus  can  provide  rotational  temperature  information  of  the  ground  state  02(X3Lg  ) 
molecules,  and  will  be  illustrated  further  in  detail  in  the  next  section. 

2.1.2  Stellar  Occulation  Spectra 

Fig.  2.5  presents  the  calculated  occultation  spectra  in  the  spectral  region  between  50000  and  56000 
cm’1  for  tangent  heights  at  60,  80  and  100  km.  Since  the  column  density  of  the  absorber,  O2,  is 


u 


Figure  2.5 
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small  for  measurements  at  high  tangent  heights,  the  overall  absorption  is  weaker  in  comparison  to 
measurements  at  lower  tangent  heights.  The  structure  of  the  Schumann-Runge  bands  can  clearly 
be  seen,  from  the  weaker  (2,0)  band  near  50000  cm1  to  the  stronger  (13,0)  band  near  56000  cm1. 
Complete  absorption  is  observed  near  56000  cm"1  for  the  measurement  at  60  km  tangent  height. 
The  envelope,  showing  a  transmission  of  -0.5  in  the  shorter  wavenumber  region,  results  from  the 
Herzberg  continuum  lying  underneath  the  Schumann-Runge  bands. 

Single  wavelength  (photometric)  occultation  measurements  carried  out  in  the  past  have  assumed  the 
knowledge  of  the  cross-section.  Because  the  absorption  cross-sections  of  the  Schumann-Runge 
bands  are  temperature  dependent,  the  accuracy  with  which  we  can  retrieve  the  O2  number  density 
profile  depends  upon  how  accurately  we  know  the  temperature  profile.  Consequently,  it  is  best  to 
carry  out  temperature  measurements  simultaneously  with  the  density  measurements.  This  can  be 
accomplished  by  carrying  out  spectrometric  occultation  measurements  in  order  to  provide 
simultaneous  spectral  information,  which  can  be  used  to  retrieve  the  temperature  profile. 

Here  we  demonstrate  the  recovery  of  the  number  density  and  temperature  profiles  by  simulating 
occultation  measurements  using  a  simple  Ebert-Fastie  Spectrometer  with  a  CCD  (Charged-Coupled 
Device)  detector.  The  intensity  signal  of  a  star,  C(zt,v),  measured  in  counts/sec  at  a  tangent  altitude 
z,  and  wavenumber  v,  can  be  expressed  by 

00 

C  (z,,  v)  =  S  J  I  (V)  T  (v‘  -  v)  dv',  (2.11) 

0 

where  S  =  AQeT0,  A  is  the  aperture  area  of  the  telescope,  Qe  is  the  quantum  efficiency  of  the 
detector,  T0  is  the  optical  transmission  of  the  system,  and  T(v'  -  v)  is  the  broadening  function  or 
transmission  function  of  the  dispersive  element  at  wavenumber  v.  Since  a  star  is  a  point  source, 
the  instrument  transmission  function  is  a  rectangular  function  with  its  spectral  width  determined  by 
how  the  dispersed  light  is  imaged  onto  the  detector  and  by  the  number  of  CCD  pixels  combined. 
The  spectral  resolution  is  thus  variable  since  the  number  of  pixels  that  can  be  combined  is  easily 
controlled. 

Fig.  2.6  illustrates  how  the  occultation  transmission  spectrum  varies  with  the  rotational  temperature 
of  the  absorber,  O2.  A  column  density  of  O2  equal  to  1022  env2  (-80  km  tangent  height)  is  used  in 
this  illustration.  Five  rotational  temperatures,  varying  from  200K  to  300K  with  an  increment  of 
25K,  are  simulated  for  two  spectrometer  instrument  widths,  50  env1  (Fig.  2.6a)  and  1000  enr1 
(Fig.  2.6b)  respectively.  The  various  Schumann-Runge  (\)',0)  bands  are  better  resolved  by  the 
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Figure  2.6 
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higher  resolution  instrument.  The  rotational  band  structure  of  each  vibrational  band  is  also  bener 
revealed,  indicating  that  the  rotational  temperature  may  be  determined  more  accurately  if  enough 
signal  were  observed.  It  is  worth  to  note  here  that  a  25K  change  in  temperature  gives  rise  to 
approximately  5  to  10  percent  change  in  the  transmission  in  certain  temperature  sensitive  spectral 
regions.  In  other  words,  an  uncertainty  of  25K  in  temperature  may  translate  into  an  error  of  ~5  to 
10  percent  in  the  deduced  column  density,  and  this  error  varies  with  the  spectral  resolution  of  the 
measured  spectra.  In  fact,  if  we  include  regions  of  the  spectra  where  the  cross-sections  are 
temperature  insensitive  and  fit  the  spectra  with  a  temperature,  the  actual  error  will  be  much  smaller. 
This  concept  is  the  basis  of  our  data  analysis  technique  used  to  retrieve  the  O2  number  density 
profile. 


2.2  Recovery  of  Density  and  Temperature  Profiles 

2.2.1.  Data  Analysis 

There  are  many  techniques  that  can  be  used  to  recover  the  number  density  and  the  temperature 
profiles  of  the  absorber  from  the  measured  stellar  occulation  spectra.  Here  we  use  a  simple 
technique  which  filters  out  the  temperature  dependence  of  the  absorption  cross-sections  first. 
Rewriting  Eq.  (2.1)  in  terms  of  t(v)  defined  in  Eq.  (2.10)  and  a  weighted  cross-section,  o(v),  one 
obtains 


00 

I  (v)  =  I„(v)  exp  Ml  n  (s)  oej  (v,  T)  ds  )  =  I«„(v)  exp  {  -  N  (02)  o  (v)  )  , 

1 


(2.12) 


where  a(v)  is 


a(v)  = 


n  (s)  aei  (v,  T)  ds’ 


I 


n  (s)  ds 


(2.13) 


The  temperature  filtering  process  can  be  accomplished  by  defining  an  effective  temperature,  T,  at 
which  the  absorption  cross-sections  best  represent  the  weighted  cross-sections.  We  can  then 
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perform  a  non-linear  least-square  fit  to  the  spectra  to  obtain  T  and  the  column  density  at  each 
tangent  ray  height.  This  is  a  good  approximation  since  the  departure  of  the  actual  o(v)  from  the 
one  given  by  T  is  very  small  and  represents  only  a  second-order  effect  to  the  total  extinction  of  the 
light.  The  error  in  the  deduced  column  density  therefore  is  negligible.  Also  if  the  atmospheric 
temperature  does  not  vary  rapidly  with  altitude,  since  most  of  the  absorption  takes  place  near  the 
tangent  point,  the  deduced  effective  temperature  becomes  very  close  to  the  one  at  the  tangent 
altitude. 


Let  Qj  denote  the  occultadon  signal  at  wavenumber  i  and  tangent  height  j,  which  is  a  function  of 
the  effective  temperature  Tj  and  O2  column  density  Nj.  The  unknowns,  Tj  and  Nj,  may  be 
calculated  from  the  measured  spectra  using  a  non-linear  least-squares  fitting  technique.  We  first 
linearize  the  problem  by  expanding  Eq.  (2.1 1)  about  empirically  determined  approximations  (T0J, 
N0j)  for  the  temperature  and  the  O2  column  density  using  only  the  first-order  terms  in  the 
expansion: 


C  (T., )  s  C.  (T  .,  N  .)  + 

‘J  J  J  U  v  °J  OJ 


~  (T  -  T  )  + 
3T  J  °J 


TTT  (N.  -  N  .). 
dN  J  °J 


(2.14) 


The  sum  of  the  squares  of  the  residuals  between  the  theoretically  computed  signal  and  the 
measured  signal  is  calculated  by 


<82>  =  V  (M.  .-C.  )2 

>J  u 


(2.15) 


where  Mjj  is  the  measured  signal  at  wavenumber  i  and  tangent  height  j.  By  minimizing  <82>,  one 
can  obtain  a  better  set  of  estimates  for  T,  and  Nj,  or  Xj  (1=1,2),  a  general  representation  of  Tj  (1=1 ) 
and  Nj  (1=2),  from  a  single  matrix  multiplication  operation: 


ITt=l 


(2.17) 
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and 


The  procedure  is  reiterated  until  the  estimated  values  converge.  The  speed  of  convergence  depends 
on  the  proximity  of  the  initial  values  for  (X^,  X0J2)  or  (Tj,  Nj)  to  the  actual  values.  The  expression 
(Eq.  2.3)  for  the  column  density  Nj  at  tangent  height  j  is  a  simple  Abel  integral  equation  which  can 
be  inverted  to  give  the  number  density  profile  of  02  (Hays  and  Roble,  1968a,b  and  Roble  and 
Hays,  1972). 


2.2.2  Statistical  Error  Analysis 

The  accuracy  in  retrieving  the  number  density  profile  from  the  occultation  intensity  data  depends 
upon  (a)  the  number  of  data  points  per  scan,  (b)  the  number  of  spectral  measurements  per  point, 
(c)  the  statistical  noise  inherent  in  the  signal,  (d)  knowledge  of  the  absorption  cross-sections,  (e) 
knowledge  of  the  satellite  position  during  occultation,  and  (f)  departures  from  spherical  symmetry. 
The  first  three  items  are  examined  in  this  chapter.  Item  (d)  has  been  discussed  in  the  previous 
section  and  is  not  included  in  the  statistical  error  analysis  performed  here.  Item  (e)  depends  on  the 
satellite  tracking  and  the  degree  of  accuracy  in  telescope  pointing.  Its  errors  will  appear  in  the 
tangent  ray  height  determination  and  will  not  be  discussed  in  this  work.  The  last  item  may  be 
important  for  species  which  change  rapidly  with  time.  In  this  analysis,  however,  a  spherical 
symmetry  is  assumed. 


Considering  that  the  occultation  data  Mjj  has  a  Poisson  random  noise,  we  can  analyze  the  accuracy 
of  the  fitted  temperature  and  the  fitted  column  density  of  02  from  the  theory  of  error  propagation. 


(2.19) 


where  AM,j  =  (Mjj)1/2. 


The  accuracy  of  the  recovered  number  density  of  O2  depends  upon  the  error  in  the  column  density 
and  the  technique  used  to  invert  the  Abel  integral.  Here,  we  just  use  a  simple  linear  matrix 
inversion  technique  to  estimate  the  error.  Writing  Abel  integral  in  a  matrix  form,  one  obtains 
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(2.20) 


N,  =  5Xv 

k 

where  nk  is  the  number  density  of  O2  at  altitude  k,  and  Wjk  is  a  geometric  weighting  function. 


t^i/2 


for  a  spherically  symmetric  atmosphere.  Here  Az  is  the  size  of  the  altitude  grid  which  the  inversion 
is  performed  on.  The  errors  in  the  inverted  O2  number  density,  Ank,  can  then  be  estimated  by 

Avi =  X  °°S AN?  ’  (2-22) 

j 

where  (Okj  is  the  element  of  the  inverted  matrix  of  W,  and  ANj  is  the  error  in  the  column  density  of 
O2  at  tangent  height  j  given  by  Eq.  (2.19). 


2.2.3  An  Example:  39_Lambda_Orion 

We  consider  an  O-Type  star,  39_Lambda_Orion,  with  a  visual  magnitude  mv  of  3.39  to  simulate 
the  occultation  sequence.  The  spectral  irradiance  of  this  star  was  measured,  varying  from 
~1.5xl02  photons/cm2/sec/A  near  40000  cm  1  to  2.5xl02  photons/cm2/sec/A  near  60000  cm4 
(Jamar,  et  al.).  A  constant  value  of  2.0xl02  photons/cm2/sec/A  is  assumed  in  our  simulation.  The 
instrument  parameters  used  are  listed  in  Table  2.1,  in  which  a  specially  coated  CCD  detector 
sensitive  to  ultraviolet  radiation,  is  adopted. 


Table  2.1 

Ebert*Fastie  Spectrometer  Instrument  Parameters 
Used  in  Stellar  Occultation  Simulation 

Ao  :  12.56  cm2  (4  cm  in  diameter) 

Tr  :  0.3 

Ck  :  0.20  (at  -2000  A) 

I.P.  :  1  second 
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Maximum  Integration  Time  (sec) 

F  i  gu  re  2.7 


During  a  stellar  occultation  sequence,  the  rate  at  which  the  absorption  path  progresses  deeper  into 
the  atmosphere  is  determined  by  the  spacecraft  orbital  velocity  and  by  the  star  location  with  respect 
to  the  orbital  plane.  In  order  to  obtain  a  given  vertical  resolution  in  the  recovered  density  profile, 
this  rate  limits  the  integration  time  allowed  for  an  occultation  measurement.  Fig.  2.7  shows  the 
time,  or  the  maximum  allowable  integration  period,  computed  for  a  star  located  on  the  orbital 
plane,  during  which  the  tangent  ray  height  changes  by  2  km  as  a  function  of  satellite  altitude.  It 
varies  from  ~0.8  second  at  60  km  tangent  altitude  to  ~  1.1  seconds  at  120  km  tangent  altitude  fora 
satellite  altitude  of  300  km,  and  it  increases  as  the  star  locates  further  away  from  the  orbital  plane. 
For  our  simulation  we  use  an  integration  period  of  1  second. 

Fig  2.8  presents  the  results  of  simulated  signal  measurements  at  tangent  ray  altitudes  of  60  km  and 
100  km  for  two  instrument  widths,  50  cm1  (Fig.  2.8a)  and  1000  cm1  (Fig.  2.8b).  Because  the 
spectral  region  where  total  absorptions  occur  moves  toward  larger  wavenumber  (shorter 
wavelength)  as  the  tangent  height  increases,  we  assume  that  the  spectrometer  steps  its  spectral 
coverage  accordingly,  150  cm1  per  every  2  km.  We  have  also  included  a  noise  of  40  counts  to 
consider  CCD  detector  read-out  error.  This  noise  value  is  insignificant  for  the  low  resolution 
measurements,  but  becomes  statistically  more  important  as  the  instrument  width  decreases. 
Although  lower  resolutions  yield  stronger  signal  levels,  they  result  in  a  smaller  number  of  spectral 
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measurements  and  less  resolution  of  the  structure  of  the  Schumann-Runge  bands.  Consequently, 
the  estimation  of  the  accuracies  of  the  recovered  density  and  temperature  requires  a  complete 
statistical  error  analysis. 


Fig.  2.9  illustrates  the  estimated  percentage  errors  in  the  recovered  density  and  temperature  profiles 
based  on  the  spectra  shown  in  Fig.  2.8.  It  shows  that  the  spectra  obtained  at  50  cm'1  instrument 
width  gives  better  determination  of  the  densities  and  temperatures,  less  than  10  percent  in  density 
and  5  percent  in  temperature.  Using  spectra  obtained  at  1000  cm1  instrument  width,  however,  we 
recover  the  density  and  temperature  not  only  less  accurately,  but  also  with  a  large  variability  with 
altitude.  This  variability  results  directly  from  the  change  in  temperature  sensitivity  as  the  spectral 
window  shfits  with  tangent  height.  Also,  as  the  tangent  height  moves  above  100  km,  the  spectral 
window  moves  towards  the  Schumann-Runge  continuum.  Since  the  Schumann-Runge  continuum 
is  essentially  temperature  independent,  this  results  in  a  gradual  unretrievable  loss  of  temperature 
information  and  an  increasing  error  in  the  deduced  temperature.  At  higher  altitudes  the  technique 
becomes  the  one  decribed  by  Hays  and  Roble  (1973),  where  only  a  single  wavelength, 
photometric  measurement  is  needed.  The  O2  column  density  can  then  be  easily  obtained  with  an 
improved  accuracy. 


2.3  Summary  and  Discussion 

A  stellar  occultation  technique  is  presented  for  the  determination  the  atmospheric  density  profile  in 
the  critical  60  to  120  km  region  from  the  O2  absorption  spectra  near  the  Herzberg  continuum, 
Schumann-Runge  bands,  and  Schumann-Runge  continum.  Because  the  absorption  cross-sections 
in  the  Schumann-Runge  bands  are  temperature  dependent,  the  data  reduction  technique  is  different 
from  the  ones  used  in  the  past.  Multi-wavelength  spectrometric  measurements  of  the  occultation 
spectra  permits  the  simultaneous  determination  of  the  atmospheric  temperature  profile. 

Using  a  simulated  CCD  Ebert-Fastie  Spectrometer,  we  have  demonstrated  the  feasibility  of  this 
technique  and  estimated  the  statistical  accuracy  in  the  recovered  density  and  temperature  profiles. 
The  errors  are  estimated  based  on  an  instrument  whose  parameters  are  described  in  Table  2.1  and  a 
source  star,  39_Lambda_Orion,  which  has  a  spectral  irradiance  of  -200  photons/cm2/sec/A  near 
the  ultraviolet  spectral  region.  We  have  shown  that  the  O2  density  profile  can  be  recovered  with  a 
statistical  accuracy  of  less  than  10  percent.  Note  that,  while  a  in  Eq.  (2.19)  is  inversely 
proportional  to  the  square  of  the  instrument  sensitivity  S  and  the  star  intensity,  Ie,  dQ/dX  and  My 
are  proportional  to  S  and  I0.  As  a  result,  the  error  in  the  recovered  column  density  is  inversely 
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proportional  to  the  square  root  of  S  and  Io.  Thus  an  occultation  experiment  performed  on  a  source 
star  other  than  39_Lambda_Orion  would  produce  different  observed  signal,  M,j,  and  results  in 
different  errors  in  the  O2  number  density  profile.  Similarily,  one  can  also  improve  the  accuracy  by 
designing  an  instrument  with  a  larger  telescope  aperture  in  order  to  increase  the  instrument 
sensitivity.  Fig.  2.10  shows  the  telescope  diameter  required  in  order  to  recover  density  with  the 
same  accuracy  demonstrated  here  for  various  star  temperatures  and  visual  magnitudes.  There  are 
approximately  500  stars  of  visual  magnitude  4.0  or  brighter  and  of  these,  about  40%  are  in  the 
spectral  classes  of  O  and  B  which  have  an  effective  temperature  greater  than  15,000K.  A  telescope 
size  similar  to  the  one  used  in  this  simulation  is  therefore  enough  to  carry  out  the  occulation 
experiment.  Of  course,  this  size  is  governed  by  the  accuracy  requirement  in  the  recovered  density 
profile  and  the  distribution  of  stars  near  the  orbital  plane  of  the  spacecraft. 

There  has  recently  been  some  interest  in  the  role  of  photodissociation  of  isotopic  oxygen  in  the 
production  of  ozone.  The  two  isotopes  of  oxygen  are  O^Ois  with  an  abundance  of  0.41%  of 
0i6016,  and  OigOp  with  -20%  of  the  Oi6Oi8  abandunce.  Although  the  isotopes  differ  slightly  in 
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mass,  they  have  quite  different  absorption  properties  from  the  ordinary  oxygen  molecules.  The 
fact  that  the  isotopes  have  twice  as  many  lines  does  not  increase  the  amount  of  absorption  because 
the  population  of  the  initial  rotational  line  also  decreases  by  about  1/2.  It  is  that  the  possibility  of 
overlapping  of  lines  and  the  variation  of  tangent  altitude  where  the  total  absorption  occurs  might 
have  some  effect  on  our  analysis.  Because  of  the  relatively  small  abundance  of  the  isotope,  we 
believe  that  the  effect  is  small.  However,  it  is  still  interesting  enough  to  be  worthy  of  further 
investigation. 

Although  we  have  only  simulated  the  occultation  process  for  the  altitude  region  between  60  and 
120  km,  this  technique  can  also  be  used  to  obtain  the  O2  density  profile  above  120  km  by 
measuring  the  absorption  in  the  Schumann-Runge  continuum  (Roble  and  Hays,  1973).  Near  60 
km  tangent  height  and  below,  however,  the  absorption  by  ozone  in  the  Hartley  bands  becomes 
important.  Ozone  has  a  peak  absorption  cross-section  of  the  order  of  10-17  cm2  near  40000  cm'1, 
giving  an  optical  depth  of  1.0  at  approximately  65  km.  Therefore,  the  technique  presented  here 
will  work  well  in  the  region  above  70  km.  Between  60  and  70  km  only  the  occultation  data  at  v  > 
45,000  cm1  a'e  useful  and  should  be  analyzed  with  extreme  care. 

There  are  many  techniques  which  can  be  used  to  retrieve  the  density  and  temperature  profiles  from 
the  occultation  data.  A  method  which  first  filters  out  the  temperature  dependence  of  the  cross- 
section  and  obtains  the  column  density  of  O2  for  each  tangent  height  was  demonstrated.  The 
systematic  error  caused  by  this  approximation  is  believed  to  be  very  small.  In  fact,  considering  the 
rate  of  data  aquisition,  we  believe  that  this  technique  is  ideal  because  of  its  simplicity  and 
efficiency. 


24 


3 .  RESONANCE  FLUORESCENCE  OF  02 


The  Atmospheric  bands  of  molecular  oxygen  originate  in  a  transition  between  the  excited  b^g* 
state  and  the  ground  X3^-  state.  The  emissions  were  first  discovered  by  Meinel  (1950)  in  the 
nightglow  spectra  shortly  after  the  discovery  of  OH  bands.  Since  then  there  has  been  an  extensive 
set  of  observational  data  under  nighttime  and  auroral  conditions.  The  (0-0)  band  of  the  02 
atmospheric  system  is  completely  absorbed  by  the  lower  atmosphere  and  can  only  be  observed 
from  altitudes  above  the  stratosphere.  As  a  result,  all  the  ground-based  measurements  have  been 
made  on  the  (0-1)  and  (1-1)  bands  (Meinel,  1951;  Chamberlain,  et  al,  1954;  Wallace  and 
Chamberlain,  1959),  and  the  (0-0)  band  emissions  were  only  observed  by  instruments  on  board 
rockets  (Packer,  1961;  Deans  et  al.,  1976;  Witt  et  al.,  1979;  McDade  et  al.,  1986)  and  the  shuttle 
(Tore  et  al.,  1985). 

Wallace  and  Hunten  (1968)  first  observed  the  dayglow  O2  (0-0)  and  (0-1)  bands  with  rocket-borne 
photometers  and  spectrometers  in  the  altitude  range  between  39-128  km.  In  their  classic  paper, 
they  presented  what  have  become  an  accepted  theory  of  the  excitation  mechanisms  for  the 
O^b'Ig*)  molecules:  resonance  fluorescence  process  in  the  region  between  70  and  100  km  and 
energy  transfer  reaction  between  O('D)  and  O2,  with  the  0(JD)  atoms  produced  in  Schumann- 
Runge  O2  dissociation  above  100  km,  and  those  produced  in  the  Hartley  dissociation  of  ozone  near 
45  km.  They  included  a  direct  contribution  from  the  photolysis  of  ozone,  which  later  was  found 
to  be  ineffective  (Gauthier  and  Snelling,  1970,  1971).  In  order  to  account  for  an  excess  emission 
near  90  km,  Wallace  and  Hunten  (1968)  also  postulated  a  small  three-body  recombination  or 
Chapman  type  chemical  source.  This  chemical  source  is  mainly  responsible  for  the  Atmospheric 
band  emissions  in  the  nightglow,  and  its  exact  mechanisms  have  been  debated  for  many  years. 
Later  studies  suggested  the  Barth  mechanism  ,  in  which  a  two-step  process  is  involved, 

O  +  O  +  M - >  O2**  +  M 


followed  by 


O2”  +  O2 - >  02(b1Eg+)  +  O2 

(Campell  and  Gray,  1973;  Greer  et  al.,  1981;  Slanger  and  Black,  1977). 

The  O2  Atmospheric  bands  are  one  of  the  brightest  emission  features  that  are  seen  in  the  visible  and 
near  infra-red  region  of  the  airglow  spectra.  The  total  band  brightness  is  strongly  coupled  to  the  02 


density  through  both  the  excitation  mechnisms  and  the  self-absorption  processes  (for  the  0-0 
band),  with  the  brightness  of  each  rotational  line  determined  by  the  atmospheric  temperature.  Here 
we  will  describe  how  the  observed  intensities  of  the  O2  Atmospheric  bands  are  used  to  retrieve  the 
O2  number  density  and  temperature  profiles  in  the  region  between  60  and  120  km. 

3.1  The  O2  Atmospheric  Bands 

3.1.1  Spectral  Intensity  of  the  Atmospheric  bands 

Fig.  3.1  shows  the  viewing  geometry  of  a  spacecraft- borne  instrument  measuring  limb 
brightnesses  of  the  O2  Atmospheric  band  emissions.  The  measured  brightness  is  an  integrated 
volume  emission  rate  along  the  slant  path  associated  with  a  particular  tangent  ray  altitude.  When 
the  absorption  is  important,  such  as  the  O2  Atmospheric  (0-0)  band,  the  measured  brightness  is 
modified  by  the  attenuation  along  the  path  between  production  region  and  the  spacecraft 

If  we  denote  tnc  total  band  volume  emission  rate  at  a  point  along  the  line  of  sight  by  ti(s),  the 


Figm  3.1 
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observed  spectral  brightness,  B(v,zt),  of  a  rotational  line  at  Vj  as  a  function  of  wavenumber  v  and 
tangent  height  z,  is  given  by 


B  (v,  Zj)  =  Jn(s)  f(s)  D(v,  s)  exp 

where  f(s)  is  the  line  strength  of  the  rotation  line  of  interest,  determined  by  the  rotational 
temperature  of  the  excited  O^b1!.*)  molecules  and  £  f  =  1.0,  D(v,s)  is  the  line  shape  at 

all  lines 

temperature  T(s),  n(s')  is  the  number  density  of  O2.  and  o(s)  is  the  self-absorption  cross-section. 
Converting  the  line-of-sight  integral  in  Eq.(3.1)  to  an  integral  in  altitude  z,  one  obtains 

B(v,  zt)  =  |  T)(z)  f(z)  D(v,  z)  W(z,  zt)  exp  {-r(v,  z,  zt)}  dz  (3.2) 

00 

+  J  T](z)  f(z)  D(v,  z)  W(z,  Zj)  exp  {-r(v,  z,  •  z,)  -  x(v,  z,  z^}  dz 

where  W(z,zt)  is  the  Jacobian  in  the  transformation  (W  =  3s/9z),  or  the  geometric  weighting 
function, 

W(z,  z,)  =  -  -  =  ,  (3.3) 

y  (Re  "*■  Z)2  '  (Re  +  Z,)2 

where  R«  is  earth  radius,  and 

00 

T(v,  z,  z,)  =  J  n(z')  o(z0  D(v,  z')  W(z',  Z()  dz',  (3.4) 

Z 

z 

t(v,  z,  zt)  =  J  n(zO  O(z0  D(v ,  z)  W(z',  zj  dz'.  (3.5) 


-Jn(s')o(s')D(v,  s^ds'  j- ds  (3.1) 


Since 
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r(v,  z,  Zj>  =  r(v,  zt,  Zj)  -  x(v,  z,  zl), 


(3.6) 


and  by  defining 

r(v,  Zj)  =  r(v,  Zj,  Zl),  (3.7) 

we  can  rewrite  Eq.  (3.2)  as 

B(v,  Zj)  =  2  J  rj(z)  f(z)  D(v,  z)  W(z,  zt)  exp  {-r(v,  Zj)]  cosh[x(v,  z,  zt)]  dz.  (3.8) 

The  above  equation  computes  the  spectral  brightness,  B(v,Zt),  of  any  given  rotational  line  when  the 
self-absorption  process  is  considered.  Since  f,  D  and  a  are  a  function  of  temperature  the  spectral 
brightness  in  this  case  not  only  depends  upon  the  band  volume  emission  rate  and  the  O2  number 
density  profiles,  but  also  on  the  atmospheric  temperature  profile.  For  emissions  that  the  self¬ 
absorption  process  is  not  important,  such  as  the  O2  Atmospheric  (0-1)  and  (1-1)  bands,  Eq.  (3.8) 
becomes  the  Abel  integral  equation  (Hays  et  al.,  1973), 

B(v,  zt)  =  J ri(z)  f(z)  D(v,  z)  W(z,  zt)  dz,  (3.9) 

and  the  spectral  brightness  is  determined  by  the  volume  emission  rate  and  the  atmospheric 
temperature  profiles. 


3.1.2  Characteristics  of  the  O2  Aunospheric  Band  Emission 
O2  Atmospheric  band  Volume  Emission  Rate 


If  we  denote  the  various  sources  of  production  of  02(b1Ig+)  by  Pres,  Po(iD)  and  Pchem  and  consider 
loss  processes  due  to  quenching  and  radiation  in  a  simple  photochemical  equilibrium  situation,  we 
can  determine  the  total  band  volume  emission  rate  as  a  function  of  altitude, 
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where  AV'v"  is  the  transition  probability  between  02(b1X,\>‘)  and  02(X3£,v»”)  states,  Ai£  is  the 
transition  probability  of  (^(b1!)  state,  and  kN2  and  k©2  are  the  quenching  rate  coefficients  of 
(^(b1!)  by  N2  and  O2  respectively. 

Fig.  3.2  shows  the  volume  emission  rate  profile  of  the  (0,0)  band  in  the  daytime  and  its  individual 
contribution  due  to  each  of  the  three  production  sources  (Bucholtz  et  al.,  1986).  Note  that  the 
dominant  sources  of  the  excited  (^(b1!)  molecules  in  the  60  to  130  km  range  are  resonance 
scattering  and  reaction  of  CH'D)  atom  with  O2  followed  by  photodissociation  of  O2.  In  a  region 
above  -90  km  where  quenching  of  (^(b'L)  molecules  is  not  important,  the  volume  emission  rate  is 
therefore  related  to  the  number  density  of  ground  state  O2.  As  the  quenching  becomes  important  in 
the  region  between  60  and  90  km,  the  production  rate  and  the  loss  rate  of  02(b1£g+)  are  both 
proportional  to  the  molecular  oxygen  number  density.  As  a  result,  the  total  band  volume  emission 
rate  remains  independent  with  altitude,  and  provides  little  information  about  the  number  density  of 
C>2.  The  capability  of  retrieving  atmospheric  density  in  this  region  lies  not  on  the  volume  emission 
rate  profile,  but  on  the  important  characteristics  of  the  (0-0)  band,  the  self  absorption  processes. 
Here  the  volume  emission  rate  profile  computed  by  Bucholtz  et  al.  (1986)  and  a  model  atmosphere 
listed  in  Table  3.1  are  used  in  our  study. 


Volume  Emission  Rate  (cm  3  sec  !) 

S-'igire  3.2 
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TABLE  3.1 


Altitudeflan) 

02(0-0)  band 

Temperature  (K)  n(0>)  n(photons/cn>Vsec) 

Contiuum  emission 
n  ( photons/cm2/sec/A) 

60.0 

252.4972 

1.49  (15) 

7323.35 

7.47  (1) 

62.0 

244.3 

1.18  (15) 

7133.66 

5.49  (1) 

64.0 

236.0 

9.26  (14) 

7089.15 

4.04  (1) 

66.0 

227.8 

7.19  (14) 

7113.96 

2.97  (1) 

68.0 

219.5 

5.54  (14) 

6934.56 

2.18  (1) 

70.0 

211.2 

4.22  (14) 

6804.00 

1.60(1) 

72.0 

208.4 

3.11  (14) 

6653.63 

1.18(1) 

74.0 

205.3 

2.28  (14) 

6553.28 

8.66  (0) 

76.0 

202.4 

1.67  (14) 

6468.28 

6.37  (0) 

78.0 

199.5 

1.22(14) 

6370.40 

4.68  (0) 

80.0 

196.6 

8.86  (13) 

6281.98 

3.44  (0) 

82.0 

193.5 

6.39  (13) 

6222.41 

2.53  (0) 

84.0 

190.1 

4.60(13) 

6278.27 

1.86  (0) 

86.0 

188.6 

3.24(13) 

6440.33 

1.36  (0) 

88.0 

186.7 

2.25  (13) 

6570.99 

1.00  (0) 

90.0 

183.8 

1.54  (13) 

6488.03 

0.73  (0) 

92.0 

186.1 

1.05  (13) 

6086.06 

0.54  (0) 

94.0 

187.9 

7.15  Cl) 

5611.43 

0.39  (0) 

96.0 

186.8 

O.'  vi2) 

5134.98 

0.29  (0) 

98.0 

189.3 

1  0(12) 

4751.40 

0.21  (0) 

100.0 

196.3 

i.92  (12) 

4424.45 

0.15  (0) 

102.0 

205.7 

1.26(12) 

4061.11 

0.11  (0) 

104.0 

216.5 

8.41  (11) 

3653.11 

8.58(-2) 

106.0 

228.4 

5.63  (11) 

3229.26 

6.30(-2) 

108.0 

240.8 

3.80(11) 

2786.08 

4.63(-2) 

110.0 

253.7 

2.59  (11) 

2309.08 

3.40(-2) 

112.0 

267.4 

1.77(11) 

1822.02 

2.50(-2) 

114.0 

283.4 

1.22  (11) 

1370.23 

1.84(-2) 

116.0 

305.1 

8.42(10) 

987.13 

1.35(-2) 

118.0 

334.7 

5.80(10) 

692.65 

9.95(-3) 

120.0 

363.7 

4.12  (10) 

490.63 

7.31  (-3) 

122.0 

391.1 

3.02  (10) 

335.93 

5.38(-3) 

124.0 

417.2 

2.27  (10) 

247.80 

3.95(-3) 

126.0 

442.1 

1.75  (10) 

186.91 

2.90(-3) 

128.0 

465.8 

1.38  (10) 

144.05 

2. 1 3(-3) 

130.0 

488.3 

1.11  (10) 

113.21 

1.57(-3) 

132.0 

509.8 

9.06  (9) 

90.53 

1.15(-3) 

134.0 

530.2 

7.49  (9) 

73.48 

8.49(-4) 

136.0 

549.7 

6.26  (9) 

60.41 

6.24(-4) 

138.0 

568.2 

5.29  (9) 

50.20 

4.58(-4) 

14J.0 

585.8 

4.50  (9) 

42.10 

3.37(-4) 
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Temperature  Effect 


The  spectral  brightness  of  any  given  rotational  line  depends  on  the  atmospheric  temperature  profile 
in  two  different  ways.  One  is  the  thermal  distribution  of  the  rotational  levels  of  both  the  emitter, 
the  excited  02(b1Zg+)  and  the  absorber,  the  ground  state  C^CX^Lg).  The  former  determines  the  line 
volume  emission  rate  profile,  while  the  latter  determines  the  line  absorption  cross-section 
(McClatchey  et  al.,  1973), 


f  =  fsYexP 


1.439  E'  (T-T.) 


TT. 


^  =  CfSYexp 


1.439  E"  (T  -  Ts) 

TT] 


(3.11) 


(3.12) 


where  fs  is  the  emission  line  strength  at  standard  temperature  Ts,  os  is  the  self-absorption  cross- 
section  at  Ts,  and  E'  and  E"  are  the  energies  in  cm1  for  the  upper  and  lower  states  of  the  rotational 
line  respectively  (Babcock  and  Herzberg,  1948). 

Fig.  3.3  presents  the  rotational  line  intensity  distributions  of  the  O2  Atmospheric  (0-1)  band  for 
temperatures  of  150  and  300*K  respectively.  The  band  origin  is  located  near  8645 A  with  the  R- 
branch  lines  appearing  on  the  shorter  wavelengths  side  and  the  P-branch  lines  appearing  on  the 
longer  wavelength  side.  Because  collisions  with  the  ambient  atmospheric  particles  are  very 
frequent  in  the  region  between  60  and  120  km,  both  the  excited  and  the  ground  states  of  O2  are  in 
rotational  thermal  equilibrium  and  have  a  rotational  temperature  equal  to  the  ambient  kinetic 
temperature.  As  the  temperature  increases,  the  lines  near  the  bandhead  and  the  bandtail, 
representing  the  transitions  from  higher  rotational  levels,  become  more  pronounced.  This  apparent 
change  in  band  structure  can  provide  us  with  the  temperature  information  of  the  atmosphere. 

If  the  self-absorption  process  is  important,  such  as  for  the  O2  Atmospheric  (0-0)  band  emission, 
the  band  strucuture  will  be  modified  due  to  variations  in  extinction  from  line  to  line.  Fig.3.4a 
shows  the  computed  intensities  of  all  the  rotational  lines  in  the  O2  Atmospheric  (0-0)  band  at  60  km 
tangent  height,  and  Fig.  3.4b  shows  the  case  when  the  self-absorption  process  is  ignored.  The 
effect  of  the  self-absorption  process  is  small  for  the  lines  near  the  bandhead  and  the  bandtail  whose 
line  strengths  are  weak,  and  is  large  for  lines  near  the  band  origin.  The  temperature  deduced  from 
the  rotational  band  structure  of  the  O2  Atmospheric  (0-0)  band  emission  therefore  is  larger  than  the 
rotational  temperature  of  the  emitters. 
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Figure  3.3 


The  other  temperature  dependence  arises  when  the  spectral  line  shapes  of  the  emission  and  the 
absorption  cross-section  are  considered.  In  the  region  above  60  km  the  pressure  broadening  effect 
is  not  very  important.  The  spectral  shapes  of  the  emission  line  and  the  self-absorption  cross- 
section  are  determined  by  the  kinetic  temperatures  of  the  excited  state  C^b1^)  and  the  ground 
state  C>2(X3Xg-)  respectively, 


D(v,  T)  =  J_ 


ln2 

~K 


i_ 


exp 


(v  -  vQ)2  ln2  ^ 


a 


D 


J 


(3.13) 


where  v0  is  the  wavenumber  at  line  center,  and  cto  is  the  doppler  width. 


v„  2kT  ln2  f 
aD  =  -l~_S— J- 


(3.14) 


Fig.  3.5  show:*  ,he  emission  line  shape  of  line  pPll  (13118.0332  cm1)  in  the  (0-0)  band  for  a 
measurement  at  a  tangent  height  of  60  km.  The  emission  line  shape  for  neglecting  the  self¬ 
absorption  process  is  also  shown.  Since  the  magnitude  of  the  self- absorption  cross-section  is  the 
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Figure  3.5 
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Brightness  (Rayleigh) 

Figure  3.6 


largest  at  the  line  center  and  decreases  toward  the  wings,  the  effect  of  the  self-absorption  process 
gives  rise  to  a  broadened  emission  line  shape  and  an  effective  temperature  larger  than  the  kinetic 
temperature  of  the  emitters. 

Oi  number  density 

One  would  also  expect  that  emission  spectral  line  shape  at  lower  tangent  heights  will  appear 
broader  than  the  ones  at  higher  tangent  heights  due  to  longer  absorption  paths  and  the  larger 
extinctions  near  the  line  center.  Consequently,  the  spectral  intensity  of  lines  in  the  (0-0)  band  not 
only  contains  the  temperature  information  of  the  emitter,  but  also  the  column  density  of  the 
absorber  O2. 

Fig.  3.4  clearly  shows  that  for  some  rotational  lines,  approximately  60-80%  of  the  intensities  are 
absorbed.  Fig.  3.6  presents  the  intensity  as  a  function  of  tangent  height  for  the  line  shown  in  Fig. 
3.5.  The  variation  of  the  amount  of  self-absorption  with  tangent  height  reflects  the  variation  of  the 
O2  column  density.  Hence,  by  measuring  the  extinction  of  the  02(0-0)  band  emission  intensity  as  a 
function  of  tangent  height,  one  can,  in  principle,  determine  the  number  density  profile  of  O2. 
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3.2  Recovery  of  Density  and  Temperature  Profiles 
3.2.1  Mathematical  formulation 

As  the  Eq.(3.8)  indicates,  the  observed  limb  brightess  B(v,Zt)  of  an  emission  line  in  the  O2 
Atmospheric  (0-0)  band  is  a  non-linear  function  of  the  volume  emission  rate,  temperature,  and  the 
O2  number  density  profiles.  Using  linear  perturbation  theory,  we  can  express  actual  n(z),  T(z), 
and  t)(z)  profiles  as  sums  of  guessed  profiles,  n0(z),T0(z)  and  ri0(z)  and  small  perturbations 
n'(z),T'(z)  and  T|'(z),  as  follows: 


n(z)  =  rto(z)  +  n'(z). 

(3.15) 

T(z)  =  T0(z)  +T(z), 

(3.16) 

Tl(z)  =  'n0(z)  +  q,(z). 

(3.17) 

We  shall  see  below  that  Eq.  (3.8)  can  be  simplified,  expressing  any  observed  limb  brightness  in 
terms  of  a  line<u  conbination  of  the  volume  emission  rate,  temperature  and  the  O2  number  density 
profiles.  The  linear  representation  of  the  limb  brightness  can  help  us  understand  the  relative 
importance  or  the  contribution  of  the  three  parameters  to  the  observed  brightness. 

Substituting  Eqs.  (3.15),  (3.16)  and  (3.17)  into  the  expressions  for  r(v,zt)  and  x(v,z,zt)  and 
neglecting  the  terms  involving  the  multiplication  of  two  perturbations  parameters,  one  obtains 

r(v,  zt)  =  r0(v,  zt)  +  rT'(v,  zt)  +  rn(v,  zt),  (3. 1 8) 

and 

T(V,  z,  zt)  =  T0(V,  z,  zt)  +  x'T(v,  z,  zt)  +  x'n(v,  z,  zt),  (3.19) 

where 


(3.20) 


(3.21) 
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and 


r„(v,  z,)  -  J  n'(z)  <70(z)  D0(v,  z)  W(z,  zt)  dz. 


x0(v,  z,  zt)  =  J  nc(z)  c0(z)  D0(v,  z)  W(z,  z,)  dz, 
z'T(v,  z,  z.)  .  Jn0(z)  r(z)  W(z,  z.)  dz, 

Zt 

z 

x'n(v,  z,  zt)  =  J  n'(z)  o0(z)  D0(v,  z)  W(z,  zt)  dz. 


The  terms  e  r  and  cosh(x)  in  Eq.  (3.8)  can  then  be  expanded  as 
e  r  =  e  r°  (1  -  r'T  -  r'n), 
cosh(x)  =  cosh(x0)  +  (t'T  +  x'n)  sinh(x0), 


and  Eq.  (3.8)  becomes 

oo 

B(v,  z()  =  B„(v,  zt)  +  j  [x'T(v,  z,  zt)  a(v,  z,  zt)  -  T T(v,  zt)  P(v,  z,  zt)]  dz 

oo 

~  j. 

+  I  [x'na(v,  z,  Zt)  -  rn  P(v,  z,  Zt)l  dz  +  JT'(Z)  tfv,  z,  z,)  dz 

oo 

+  |  rj'(z)  5(v,  z,  zt)  dz, 

A 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


(3.26) 

(3.27) 


(3.28) 


where 


B0(v,  zt)  =  J  2n0(z)  fo(z)  D0(v,  z)  W(z,  z()  exp  {-ro(v,  zt)}  cosh  [x0(v,  z,  zt)]  dz, 


(3.29) 


and 
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a(v,  z,  Zj)  =  2rj0(z)  fQ(z)  D0(v,  z)  W(z,  zj  exp  {-ro(v,  z,)}  sinh  [x0(v,  z,  z,)],  (3.30) 

P (v,  z,  Zj)  =2r[0(z)  fe(z)  D0(v,  z)  W(z,  Zj)  exp  {-ro(v,  zt))  cosh  [t0(v,  z,  z,)],  (3.31) 


7  (v,  z,  Zj)  =  p(v,  z,  zt) 


1  3Do  +.  1 


D0(v,  z)  3T  f0(z)  9T 
8(v,  z,  z,)  =  P(v,  z,  zt)  /  T)0(z). 


(3.32) 

(3.33) 


Here,  the  derivatives  of  f(T)  and  D(v,T)  with  respect  to  temperature  T  can  be  easily  obtained  from 
Eq.  (3.11)  and  (3.13), 


IiL=1  T 

f  ar  l  T 

_L  ao  _  j_{ 
D  9T  *  2T  | 


-1+1-439  E' 


r2 


2  ln2  (v  -  v0)2  1 


a 


J 


(3.34) 


(3.35) 


Expressing  the  integration  in  Eq.  (3.28)  for  a  given  tangent  height  measurement  i  in  a  form  of 
matrix  summation,  one  obtains 


(  '  \ 
\ 

— 

V  *  / 

k=I 

k=l 

V  *  / 

(3.36) 


where 


B, 


fc=i 


Z*=IV 


£vS> 

L  >* 

£  v  1  p„ 


°kDk 


+  4’ 


(3.37) 

(3.38) 

(3.39) 

(3.40) 
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and 


1  3(oD)  _  1  .  1.439  E"  1 

oD  3T  7-^27 


2  In2(v  -  Vp) 

a  i 


- 1 


(3.41) 


Note  that  Bj’(v)  is  wavenumber  dependent,  and  for  an  instrument  with  transmission  function  I(v- 
v0)  and  sensitivity  S  the  perturbed  signal  in  counts/sec  at  wavenumber  v0  can  be  calculated  as 


Bi<vo)  =  BJV.)  +  £ [\  ) X*  +  £  V Y*  +  £ 

le=l  v  K  '  lc=l  lc=l  v  ' 


"ilc 


where 


(3.42) 


» 

\  ®io  (V)  = 


I(v  -  v0)  dv. 


and 


X^v.)  =  S  J  X*  I(v  -  v0)  dv, 

O 


(3.43) 


(3.44) 


Y&  (v0)  =  S  J  Yj^  I(v  -  vQ)  dv, 

O 


(3.45) 


\  (v0)  =  S  J  I(v  -  v0)  dv. 

O 


(3.46) 


3.2.2  Data  Analysis 

Consider  B,(v0)  as  the  measured  signal  at  wavenumber  vQ  and  Bj0(v0)  as  the  theoretically  estimated 
signal  based  upon  a  guessed  volume  emission  rate  profile  t|0,  a  tempeature  profile  T0(z),  and  an 
number  density  profile  nQ(z).  Eq.  (3.42)  therefore  determines  B;(v0)-  B;0(v0),  the  deviation  of  the 
observed  signal  from  the  estimated  signal  based  upon  the  three  guessed  profiles. 

The  X,Y,  and  Z  in  Eq.  (3.42),  in  general,  are  referred  to  as  the  weighting  functions.  They 
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respresent  the  weightings  of  each  parameters  from  various  altitudes,  in  this  case,  the  n'/n0,  T,  and 
tj'AIo  respectively,  to  the  difference  between  the  observed  and  the  computed  signals.  If  X  is  the 
largest  among  the  three,  the  observed  deviation  is  most  sensitive  to  the  guessed  n(z)  profile, 
making  the  O2  number  density  the  best  parameter  to  be  recovered. 

There  are  many  techniques  that  may  be  used  to  recover  ri(z),  T(z)  and  n(z)  from  the  measurements 
of  Bj(v).  One  is  the  least-square  fitting  technique,  in  which  iteration  is  performed  to  find  the  best 
tj(z),T(z)  and  n(z)  profiles  that  give  the  minimum  deviation  B{(n0)-  Bi0(n0)  summed  over  all 
spectral  measurements.  Since  there  are  three  sets  of  profiles  needed  to  be  recovered,  multi¬ 
wavelength  spectroscopic  measurements  (at  least  three  spectral  region)  of  the  O2  Atmospheric  (0-0) 
band  are  required.  The  data  analysis  task  involved  in  this  technique  is  very  inefficient  and  requires 
inversions  of  matrices  in  large  dimensions.  The  large  number  of  parameters  involved  in  the  fitting 
process  not  only  makes  the  iteration  convergence  difficult,  but  also  produces  large  errors  in  the 
recovered  profiles. 

Another  alternative  is  to  measure  temperature  and  volume  emission  rate  profiles  simultaneouly  with 
the  O2  Atmospheric  (0-0)  band  emission.  This  can  be  accomplished  by  carrying  out  spectroscopic 
measurements  of  the  (0-1)  band  emission.  The  volume  emission  rate  of  the  O2  Atmospheric  (0-1) 
band  is  linearly  related  to  that  of  the  (0-0)  band  with  a  proportional  constant  determined  by  the 
Franck-Condon  factor.  Because  the  (0-1)  band  emission  is  optical  thin,  the  measured  line-of-sight 
brightness  is  a  simple  Abel  integral  and  its  volume  emission  rate  profile  can  be  retrieved  efficiently 
(Hays  et  al.,  1973).  In  addition,  the  temperature  can  also  be  obtained  from  the  spectroscopic 
measurements  of  the  rotational  structure  of  the  (0- 1 )  emission  band. 


(0-1)  band:  Temperature  and  volume  emission  rate  measurements 

If  the  (0-1)  band  emission  is  measured  by  a  low  resolution  spectrometer  whose  instrument  width  is 
much  wider  than  the  doppler  width  of  the  emission  line,  the  observed  signal  at  wavenumber  v0  is 

00 

B(v0,  zt)  =  £sf  ^(z)  fj(z)  W(z’  Z‘}  I(VJ '  Vo)  dz’  (347) 

**  l 

where  v,  is  the  wavenumber  of  the  rotational  line  j,  and  I(vj-v0)  and  S  are  the  instrument 
transmission  function  and  the  sensivity  respectively. 


40 


Fig.  3.7  shows  the  computed  signals  of  the  O2  (0-1)  band  emission  using  a  12-channel  Image 
Plane  Detecter  Spectrophotometer  (IPDS)  (Appendix  I)  for  two  tangent  height  measurements:  60 
and  80  km.  A  volume  emission  rate  profile  derived  from  the  one  for  the  (0-0)  band  (Bucholtz  et 
al.,  1986)  and  a  continuum  scattering  background  which  decreases  exponentially  with  altitude  have 
been  assumed.  The  instrument  parameters  of  the  IPD  Spectrophotometer  listed  in  Table  3.2  and  a 
field-of-view  corresponding  to  an  image  of  2km  in  the  vertical  and  300  km  in  the  horizontal 
directions  of  the  Earth's  limb  are  used  in  the  simulation. 


TABLE  3.2 

IPD  Spectrophotometer  Instrument  Parameters 
Used  in  O2  Resonance  Fluorescence  Simulation 

Ao  :  45.603  cm2  (3  inch  in  diameter) 

W  :  1.5  x  10"4  ster  (-2.0  km  x  300  km  at  the  limb) 

Tr  :  0.3 

Qe  :  4%  (at  7620  A)  and  2%  at  (8650  A) 

I.P.  1  second 


The  volume  emission  rate  at  wavenumber  vQ  and  altitude  z,  q(vD,z),  defined  as 

q(v0,  z>  =  s  X  Tl(z)  fj(z)  I(Vj  '  Vo)’  {3  48) 


can  then  be  obtained  by  using  the  Abel  inversion  algorithm  (Hays  et  al.,  1973).  After 
measurements  at  all  the  wavenumbers  have  been  inverted,  we  can  get  the  spectral  information  or 
the  rotational  band  structure  of  the  (0-1)  band  emission  rate. 

Fig.  3.8  gives  the  volume  emission  rates  recovered  at  the  12  spectral  regions  of  the  IPDS  at  60  and 
80  km  altitudes,  revealing  clearly  the  rotational  emission  band  structures.  By  using  a  nonlinear 
least-square  fitting  method  (Killeen  and  Hays,  1984;  Appendix  I),  one  can  then  obtain  the  total 
band  volume  emission  rate  and  the  rotational  temperature  of  the  02(b1Lg+)  molecules.  In  order  to 
conserve  computation,  wc  first  linearize  the  problem  by  expanding  Eq.  (3.48)  about  a  set  of 
approximations  or  initial  guesses  (q,  T,  qc)  and  keep  only  the  first-order  terms  in  the  expansions: 
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Figure  3.7 
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^kn^k,Tk,Tlck^  S  ^kn^  ko’^ko’^  ck^  ^  ^k' ^  ko^+  ^  ^k'^ko^  ^ck'^ck^  (3.49) 

where  Tjiu,  is  the  (0-1)  band  volume  emission  rate,  at  wavenumber  channel  n  and  Tic*  is  the 
continuum  scattering  rate  at  altitude  k  respectively. 


If  the  measured  volume  emission  rate  at  the  n-th  channel  obtained  from  Eq.  (3.49)  is  expressed  by 
Qm,  the  sum  of  the  square  of  the  residuals  between  the  theoretically  computed  signal  and  the 
measured  signal  Ct,  is  calculated  by 


<5k2>'Z((i»-V 


(3.50) 


2  1 
By  minimizing  <S£>,  one  can  obtain  a  set  of  better  estimated  T|k>Tk  and  rick,  or  Pfc  (1=1, 2, 3  for 

q,T,qc  respect, .  ely),  a  general  representation  of  the  newly  determined  parameters  from  a  single 

matrix  multiplication  operation: 


^ 


(3.51) 


where  Mkc  is  an  inverse  matrix  of  L^, 
t1  _  f’  ^kn  ^kn 

^  tf  3P?  <  ’ 


(3.52) 


'•  1 da 


tt 


(3.53) 


Iterations  are  performed  in  order  to  obtain  convergence  of  the  results.  The  speed  of  convergence 
depends  upon  the  proximity  of  the  initial  guessed  values  to  the  actual  values. 
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(0-0)  Atmospheric  band:  Density  Measurement 


After  the  (0-1)  band  volume  emission  rate  and  the  atmospheric  temperature  profiles  are  determined, 
we  can  calculate  from  Eq.  (3.42),  B;(v0)-  BiO(v0),  the  deviation  of  the  observed  signal  at  tangent 
height  i  and  wavenumber  vQ  for  the  (0-0)  band  emission  from  the  estimated  signal  based  upon  a 
guessed  O2  number  density  profile  n0(z). 


B'1(V0)  =  Bi(v0)-Bio(v0) 


■m 


(3.54) 


We  can  solve  the  above  linear  equation  for  n’(z)  to  get  a  better  estimate  of  n(z),  n(z)  =  n<,(z)  +  n’(z) 
bv  minimi'nng  the  square  of  the  deviation  of  B;  (v0)-  B;0(v0)  summing  over  all  spectral  and  tangent 
heights.  The  n(z)  profile  is  then  used  as  a  new  guess  to  Eq.(3.54),  and  the  process  is  repeated 
until  n’(z)  becomes  very  small  compared  to  nc(z).  This  technique  is  identical  to  the  one  used  for 
the  (0-1)  band  data  analysis  and  gives  one  best  O2  number  density  profile  based  upon  the 
measurements  at  all  of  the  tangent  heights  and  the  spectral  regions. 


Because  the  O2  number  density  information  can  only  be  recovered  if  there  is  significant  degree  of 
extinction  in  the  measured  intensity  of  the  O2  Atmospheric  (0-0)  band  emission,  the  capability  of 
this  technique  is  therefore  limited  by  the  magnitude  of  the  absorption  cross-section.  In  other 
words,  the  O2  number  density  profile  can  only  be  directly  measured  up  to  an  altitude  of 
approximately  90  km.  Since  the  temperature  profile  is  obtained  in  the  whole  region  between  60 
and  120  km,  the  O2  density  from  90  km  to  100  km  at  the  turbopause  can  be  derived  from  the 
barostatic  equation  (Banks  and  Kockarts,  1973), 


(3.55) 


under  perfect  mixing  condition,  and  above  100  km. 
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under  diffusive  equilibrium  condition,  where  m  and  mo2  are  the  mean  molecular  weights  of  the 
atmosphere  and  O2  respectively. 
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Fig.  3.9  shows  the  simulated  signals  that  the  EPDS  measures  for  a  one  second  integration  time  at  a 
tangent  height  of  60  km.  The  unattenuated  signals,  derived  from  the  (0-1)  volume  emission  rate 
profile,  are  also  shown.  Although  it  is  not  a  practical  approach,  we  can  still,  in  principle,  recover 
the  n(z)  profile  directly  from  Eq.  (3.54)  using  measurements  from  a  single  channel.  However, 
because  the  optical  depth  differs  for  each  channel,  the  accuracy  in  the  recovered  density  profile  is 
different. 


3.3  Statistical  Accuracy  Analysis 


As  indicated  in  Fig.  (3.9),  the  difference  in  the  observed  and  computed  signal,  Bi  -  Bi0,  depends 
upon  not  only  the  O2  number  density  profile  n/no,  but  also  the  temperatue  T'  and  the  volume 
emission  rate  profiles  T|'/r|0.  The  uncertainties  in  the  temperature  and  the  volume  emission  rate 
determined  from  the  (0-1)  band  emission  would  give  rise  to  error  in  the  recovered  O2  number 
density. 


Statistical  Errors  in  Temperature  and  Volume  Emission  rate 

The  statistical  errors  AP^  in  the  volume  emission  rate  (P^)temperature  (P^)  and  obtained  from  the 
(0-1)  band  measurements  can  be  estimated  from  Eq.  (3.51)  using  theory  of  error  propagation. 
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where  AC^,  =  (W^,  assuming  that  measured  signal  Ciu,  has  a  Poisson  random  noise.  Fig.  3.10 
presents  the  statistical  accuracies  in  the  temperature  and  total  (0-1)  b:  nd  volume  emission  rate 
profiles  obtained  based  upon  the  simulated  measurement'’  given  in  Fig.  3.8. 

Statiscal  Errors  in  O2  number  density 

The  capability  of  recovering  O2  number  density  profile  from  the  (0-0)  band  emission  depends  upon 
the  relative  magnitude  of  its  weighting  function  (X)  in  Eq.  (3.42)  compared  to  the  ones  of 
temperature  (Y)  and  the  volume  emission  rate  (Z).  For  example,  consider  a  case  when  the 
weighting  function  for  the  volume  emission  rate  r|’/rj0  is  the  largest  among  the  three.  A  small  error 
in  the  profile  of  q(z)  would  produce  a  large  error  in  the  computed  signal  Bio,  which  in  turn  is 
propogated  to  the  recovered  O2  number  density  profile. 
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Figure  3.10 


Fig.  3. 1 1  shows  the  weighting  functions  X,  Y,  and  Z  as  a  function  of  altitude  lor  measurements  at 
60  km  tangent  height  for  both  a  weak  (channel  3)  and  a  strong  (channel  9,  absorption  case 
respectively.  The  maximum  weighting  functions,  especially  for  the  weak  absorption  case,  are 
found  at  the  tangent  altitude,  resulting  directly  from  the  geometric  enhancement  effect.  For  the 
stronger  absorption  case,  most  of  the  emission  originating  from  the  tangent  altitude  region  is 
absorbed  and  cannot  be  observed.  As  a  result,  the  altitude  that  the  O;  number  density  can  be 
recovered  is  higher  and  the  weighting  function  near  the  tangent  altitude  are  reduced. 


Expressing  n'/n0  in  terms  of  B’,T',  and  r|'/r|0. 
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we  can  estimate  the  error  in  the  recovered  n'/n0  using  theory  of  error  propogation, 
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where  AB;  is  the  error  in  the  measured  signal,  Alnk'/njc),  AOV/rik)  are  the  relative  error  in  the  O2 
number  density  and  the  (0-0)  band  volume  emission  rate,  and  AT*’  is  the  error  in  the  temperature 
respectively. 

Fig.  3.12  shows  the  accuracies  of  the  retrieved  O2  number  density  profile  based  on  the  (0-0)  band 
12-channel  IPDS  intensity  signals  shown  in  Fig.  3.9a.  The  accuracies  of  the  temperature  and 
volume  emission  rate  profiles  recovered  from  the  (0-1)  band  measurements  have  been  included. 
The  measured  signals  at  channel  12  (-7570A)  have  been  used  as  the  contribution  of  the  scattering 
continuum  background.  Among  the  other  eleven  spectral  regions  measured,  four  are  able  to 
provide  the  O2  number  density  profile  with  an  accuracy  less  than  10%  below  90  km.  Because  both 
the  emission  signals  and  the  self-absorption  cross-sections  are  small  in  the  spectral  regions  near  the 
bandhead  (channel  11)  and  the  bandtail  (channel  1),  the  accuracies  in  the  recovered  O2  number 
density  are  very  poor.  As  the  altitude  increases  above  90  km,  the  self-absorption  process  becomes 
less  important,  resulting  in  a  rapid  deterioration  in  the  quality  of  the  recovered  O2  number  density 
profile. 


Percentage  Error  in  the  Recovered  Density 


Figure  3.12 
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Percentage  Error  in  the  Recovered  Density 

Figure  3.13 

Fig.  3.13  presents  the  accuracy  of  the  recovered  O2  number  density  profile  when  the  measurements 
from  all  the  eleven  channels  are  combined.  It  demonstrates  that  using  a  low  resolution 
spectrometer,  such  as  the  IPDS  used  here,  one  can  retrieve  O2  number  density  profile  below  90  km 
with  an  accuracy  of  less  than  5%  for  one  second  integration  time  at  each  tangent  height.  If  we  use 
the  barostatic  equation  to  constrain  the  density  profile  to  follow  the  temperature  profile,  the 
accuracy  is  improved,  especially  in  the  region  above  90  km  where  the  accuracy  from  the  direct 
measurements  is  poor. 

3.4.  Summary  and  Discussion 

A  technique  used  to  determine  the  atmospheric  temperature  and  density  profiles  in  the  critical  60  to 
120  km  region  from  the  O2  (0-0)  and  (0-1)  Atmospheric  band  emissions  is  presented.  The 
technique  requires  spectrocopic  measurements  of  the  (0-1)  band  emission  to  obtain  the  volume 
emission  rate  profile  of  the  Atmospheric  band  and  the  atmospheric  temperature  profile.  The  (0- 1 ) 
band  emission  is  optically  thin,  and  it  relates  to  the  (0-0)  band  emission  by  the  Franck-Condon 
factor.  The  O2  number  density  profile  is  then  obtained  from  the  amount  of  self-absorption 
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determined  from  the  measured  (0-0)  band  emission  brightness  and  the  unabsorbed  emission 
brightness  computed  from  the  (0-1)  volume  emission  rate  profile. 


Here  we  have  discussed  a  technique  which  is  based  upon  daytime  measurements  of  the  O2 
Atmospheric  band  emissions.  In  principle,  this  technique  should  also  apply  to  nighttime 
measurements,  although  the  emission  intensities  are  much  weaker,  approximately  1/20  of  the 
intensities  in  the  daytime.  The  O2  Atmospheric  bands  at  night  are  produced  by  the  chemical  source 
mainly  and  have  a  emission  layer  peaked  at  approximately  95  km.  Using  the  data  analysis 
technique  discussed  in  the  text,  the  temperature,  O2  density  and  the  volume  emission  rate  profiles 
can  only  be  obtained  in  a  region  between  80  and  110  km.  The  O2  density  below  80  km  therefore 
could  only  be  obtained  using  the  least-square  fitting  techniques  to  get  the  temperarnre  and  O2 
density  profile  simultaneously  from  the  amount  of  self- absorptions. 

Measurements  at  low  tangent  heights  also  include  the  scattered  sunlight  due  to  Rayleigh  and 
aerosol  scattering  processes.  The  contribution  due  to  these  two  processes  are  highly  variable, 
depending  upc ■>  the  atmospheric  and  aerosol  density  profiles,  the  scattering  angle,  and  the  effective 
surface  albedo.  Among  the  factors  which  affect  the  continuum  brightness,  the  scattering  angle  is 
the  only  one  that  we  can  control  during  measurement  sequences.  The  scattered  light,  in  general,  is 
distributed  equally  in  the  forward  and  backward  directions  for  the  Rayleigh  scattering  process  and 
concentrated  mostly  in  the  forward  direction  for  the  aerosol  scattering  process.  In  order  to 
minimize  the  impact  of  the  scattered  sunlight  on  the  accuracies  of  our  measurements,  the  instrument 
should  have  a  line-of-sight  which  points  at  90  degree  away  from  the  Sun,  01  has  a  solar  azimuthal 
angle  of  90  degree.  The  continuum  brightnesses  included  in  our  study  here  are  obtained  using  a 
single  scattering  model  developed  at  the  University  of  Michigan's  Space  Physics  Research 
Laboratory.  A  solar  azimuthal  angle  of  90  degree  is  assumed. 

We  use  a  low  resolution  12-channel  IPD  Spectrophotometer  described  in  Appendix  I  to  simulate 
the  measurements  of  the  O2  Atmospheric  band  emissions  from  space.  We  demonstrate  in  this 
report  that  for  one  second  integration  time  in  each  tangent  height  measurement,  one  can  recover  the 
temperature  and  the  O2  number  density  profiles  with  statistical  accuracies  of  less  than  2%  at  a 
vertical  and  horizontal  resolution  of  2  km  and  300  km  respectively.  Because  the  quantum 
efficiencies  of  the  IPD  are  small  at  the  spectral  regions  near  7600  and  8600A,  a  3.0-inch  diameter 
telescope  is  adopted  in  order  to  achieve  the  required  signal  levels.  As  the  ring-type  Charge- 
Coupled  Device  (CCD)  detector  becomes  available,  the  quantum  efficiency  can  be  increased  by  at 
least  a  factor  of  10.  The  size  of  the  telescope  in  this  case  can  either  be  reduced  or  be  maintained  to 
further  improve  the  accuracies  of  the  recovered  density  and  temperature.  Most  important  of  all, 
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this  would  make  the  nighttime  measurements  of  the  temperature  and  density  profiles  more  accurate. 


The  (0- 1 )  band  emission  was  used  in  our  example  to  determine  the  atmospheric  temperature  profile 
between  60  and  120  km.  Its  emission  brightness,  however,  becomes  weaker  as  the  tangent  height 
increases  above  100  km,  resulting  in  a  rapid  deterioration  in  the  accuracy  of  the  recovered 
temperature  profile  (Fig.  3.7).  In  order  to  improve  the  accuracy,  one  can  replace  the  (0-1)  band 
emission  by  the  stronger  (0-0)  emission  which  is  optically  thin  above  100  km. 

The  02(b1Lg+ )  excited  molecules  are  produced  mainly  by  the  energy  transfer  of  O0D)  atoms  with 
O2  at  altitudes  above  100  km,  following  the  photodissociation  of  O2  by  the  absorption  in  the 
Schumann-Runge  continuum.  By  knowing  the  temperature  and  the  O2  number  density  profiles  in 
the  region,  one  can  quantitavely  examine  the  role  which  this  important  photodissociation  process 
plays  in  the  Earth's  energy  budget. 
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4.  RAYLEIGH  SCATTERING  TECHNIQUE 


The  radiation  that  emerges  from  a  planetary  atmosphere  is  determined  by  a  number  of  factors 
which  depend  on  the  character  of  the  atmosphere,  the  nature  of  the  planetary  surface,  and  the  solar 
spectrum.  The  atmospheric  panicles  responsible  for  scattering  cover  a  range  of  sizes  from  gas 
molecules  (— lO8  cm)  to  large  raindrops  and  hail  particles  (~  cm).  The  relative  intensity  of  the 
scattering  light  depends  strongly  on  the  ratio  of  the  particle  size  to  the  wavelength  of  the  incident 
radiation.  When  this  ratio  is  small,  the  scattered  light  is  distributed  equally  into  the  forward  and 
backward  direction  (Rayleigh  scattering).  When  the  particles  are  large,  an  increasing  portion  of  the 
light  is  concentrated  in  the  forward  direction.  In  this  case  (Mie  scattering),  the  distribution  of 
scattering  light  intensity  with  scattering  angle  becomes  very  complex. 

The  major  gases  of  the  atmosphere  scatter  light  due  to  the  processes  of  Rayleigh  scattering.  In  this 
case,  the  scattering  coefficient  is  inversely  proportional  to  the  fourth  power  of  the  wavelength  of 
the  incident  radiation.  During  the  daytime  period  the  blue  sky  results  from  this  preferential 
scattering  at  shorter  wavelengths.  If  the  absorption  processes  are  not  important,  the  intensity  of  the 
scattering  light  contains  the  information  of  the  number  density  of  the  atmospheric  gases.  Thus 
Rayleigh  scattering  can  be  used  to  determine  the  density  of  the  atmosphere  at  high  altitudes.  This 
method  of  observation  is  illustrated  in  Fig  4.1  and  have  been  used  or,  the  Solar  Mesosphere 
Explorer  Satellite  to  determine  the  density  and  temperature  in  the  upper  stratosphere  (Rusch,  et  al., 
1983). 

However,  questions  regarding  this  method  of  measuring  the  density  arise  due  to  the  presence  of 
aerosols  and  Polar  Mesospheric  Clouds  (Thomas  and  McKay,  1985).  Rayleigh  theory  of 
scattering  cannot  be  applied  to  these  particles,  because  the  size  of  the  aerosols  are  of  the  same  order 
of  magnitude  or  larger  than  the  wavelength  of  the  incident  light.  The  basic  theory  for  the  study  of 
scattering  light  by  aerosols  was  presented  by  Mie  (1908).  This  theory  generally  assumes  an 
ensemble  of  identical  spherical  particles  and  requires  considerable  computing  resources  for  its 
solution.  Simplified  formula  have  often  been  used  (Van  de  Hulst,  1957). 

The  capability  of  measuring  the  density  of  the  atmospheric  gases  using  Rayleigh  scattering 
technique  depends  upon  the  relative  magnitude  of  the  aerosol  scattering,  since  it  causes  scattering 
that  is  difficult  to  distinguish  from  density  variations  in  the  basic  atmosphere.  We  will  address  this 
issue  in  the  next  section. 
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Figure  4. 1 


Another  question  regarding  this  technique  arises  due  to  the  variability  of  the  Earth's  surface  albedo. 
As  illustrated  in  Fig.  4.1,  the  atmospheric  particles  along  the  line-of-sight  scatters  light  which 
comes  from  two  distinct  sources.  One  comes  directly  from  the  sun  and  in  this  case,  the  emission 
rate  is  only  a  function  of  the  solar  flux,  number  density  of  the  scattering  volume,  and  the  scattering 
phase  angle.  The  other  comes  from  the  indirect  solar  radiation  which  is  reflected  by  the  Earth's 
surface.  The  relative  importance  of  this  source  depends  on  the  surface  albedo  and  the  tangent  ray 
height  of  the  measurement.  As  the  surface  albedo  varies,  the  observed  line-of-sight  intensity 
changes.  As  a  result,  this  techique  might  interpret  this  perturbation  as  density  variations  in  the 
atmosphere. 

In  this  report,  we  simulate  the  Rayleigh  and  aerosol  scattering  in  the  50-120  km  altitude  region 
using  a  single  scattering  program  generated  here  at  the  University  of  Michigan's  Space  Physics 
Research  Laboratory.  This  single  scattering  model  assumes  that  an  incoming  photon  is  scattered  at 
most  once  by  atmospheric  major  gases  (Rayleigh)  and  aerosols  (Mie).  The  effect  of  multiple 
scattering  by  underlying  layers  and  reflection  at  the  surface  are  simulated  by  numerous  effective 
albedos.  We  will  show  how  we  can  resolve  the  problem  due  to  the  existence  of  aerosol  by  using  a 
multiplicity  of  wavelengths.  We  will  also  demonstrate  the  difficulty  caused  by  the  uncertainty  in 
the  surface  albedo. 
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4.1  Scattering  Intensity  Simulation 


Effect  of  Aerosols 

Fig.  4.2  shows  the  calculated  scattering  brightness  as  a  function  of  tangent  height  due  to  Rayleigh 
only,  aerosol  only,  and  Rayleigh  plus  aerosol  scattering  for  wavelengths  at  4000A,  5000A,  6000A 
and  8000A  respectively.  The  results  shown  here  are  simulated  for  a  spacecraft  located  at  an 
altitude  of  450  km,  solar  zenith  angle  of  45  degree  and  a  surface  albedo  of  0.2. 

The  presence  of  aerosols  in  the  mesosphere  significantly  affects  the  scattered  brightness  as  a 
function  of  wavelength.  In  general  Rayleigh  scattering  takes  place  more  efficiently  at  shorter 
wavelength,  varying  with  1  A4.  The  aerosol  scattering,  however,  goes  as  lA.n,  where  n<4,  and 
increases  in  importance  at  longer  wavelength.  This  is  readily  apparent  in  Fig.  4.2.  At  4000A,  the 
contribution  to  the  brightness  profile  by  aerosol  scattering  is  an  order  of  magnitude  less  than  the 
contribution  by  Rayleigh  scattering.  The  measured  intensities  in  this  short  wavelength  spectral 
region  are  therefore  mostly  produced  by  the  Rayleigh  scattering  processes  and  are  best  suited  for 
the  determination  of  atmospheric  densities. 

For  our  calculations  we  have  used  an  aerosol  density  profile  contained  in  the  LOWTRAN  5  model. 
The  LOWTRAN  5  model  allows  the  specification  of  aerosol  levels  (i.e.  background,  low, 
moderate  and  high  volcanic  activity,  etc.)  in  the  lower  stratosphere  and  above  which  it  allows  the 
aerosol  density  profile  to  drop  off  exponentially  with  height.  The  question  arises  as  how  well  we 
know  the  aerosol  concentration  in  the  atmosphere. 

Few  observations  and  studies  of  mesospheric  aerosol  concentrations  and  distributions  have  been 
made  since  a  spate  of  observations  in  the  early  1970's.  Lidar  shows  promise  as  a  technique 
capable  of  routine  measurements  of  mesospheric  aerosol  density.  Radiance  measurements  from 
the  SME  satellite  are  a  potential  database  for  mesospheric  serosol  studies.  Aerosol  layers  from 
volcanic  and  unidentified  sources  in  the  middle  stratosphere  have  been  readily  observed  by  the 
SME  satellite  (Clancy,  1986).  A  significant  enhancement  in  aerosol  scattering  over  Rayleigh 
scattering  occurred  after  the  eruption  of  EL  Chichon  in  1982  (Naudet  and  Thomas,  1987).  A 
preliminary  SME  data  analysis  indicates  that  for  equatorial  regions  an  enhancement  in  brightness  of 
10-15%  is  due  to  mesospheric  aerosols,  and  that,  if  concentrations  of  aerosols  exists  in  the 
mesosphere  at  midlatitudes,  the  aerosols  are  well-mixed.  Clouds  in  the  summer  polar  mesosphere 
display  an  enhancement  in  brightness  in  SME  data  that  is  one  or  two  orders  magnitude  greater  than 
the  brightness  due  to  Rayleigh  or  aerosol  scattering  (Thomas  and  McKay,  1986). 
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Figure  4 . 3 


While  a  climatology  and  height  distribution  of  mesospheric  aerosols  is  not  presently  available, 
observations  suggest  that  aerosols  are  a  minor  constituent  in  the  mesosphere.  Although  the 
presence  of  a  larger  than  expected  concentrations  of  aerosols  in  the  mesosphere  would  significantly 
and  negatively  impact  on  the  Rayleigh  scattering  measurement  technique  of  atmospheric  density, 
the  judicious  selection  of  wavelength  and  strategy  can  still  minimize  the  effects  of  aerosols. 

Surface  Albedo  Effect 

Fig.  4.3  shows  the  same  results  for  Fig.  4.2  except  that  the  surface  albedo  is  increased  to  0.6.  As 
the  surface  albedo  increases,  the  contribution  to  the  scattered  brightness  from  the  ground  reflection 
also  increases.  For  example,  a  brighntness  of  4.2xl05  Rayleigh/A  at  4000A  is  obtained  for  a 
measurement  at  60  km  tangent  height  and  an  effective  surface  albedo  of  0.2.  As  the  surface  albedo 
is  increased  to  0.6  while  keeping  the  atmospheric  condition  and  the  observational  geometry  the 
same,  an  increase  of  approximately  50%  to  6.0xl05  Rayleigh/ A  in  the  scattered  brightness  is 
obtained.  For  measurements  at  120  km  tangent  height  this  brightness  also  changes  from  3.3X101 
to  4.9X101  Ravleigh/A.  The  effective  surface  albedo  is  highly  variable  due  to  the  types  of  surface 
(continent  or  ocean,  forest  or  ice  covered,  etc)  and  the  amount  of  cloud  coverage.  The  atmospheric 
density  profile  determined  from  the  scattered  intensity  therefore  is  very  uncertain. 


4.2  Summary  and  Discussion 

The  Rayleigh  scattering  technique  does  not  appear  to  be  a  viable  remote  sensing  technique  for  the 
determination  of  atmospheric  density.  This  is  due  to  the  lack  of  distinct  differentiation  between  the 
contribution  to  the  observed  brightness  from  the  Rayleigh  scattering  and  the  contributions  from  the 
aerosol  scattering  and  ground  reflection. 
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5.  MERIT  EVALUATION 


Three  techniques  which  have  the  potential  to  measure  the  atmospheric  density  in  the  altitude  region 
between  60  and  120  km  have  been  examined  Each  of  these  techniques  has  inherent  advantages  as 
well  as  weakness.  In  this  section  we  compare  the  relative  merits  of  these  related  techniques. 

5.1  Feasibility 

Rayleigh  Scattering  Technique 

•  not  feasible 

-  due  to  its  lack  of  distinct  differentiation  between  the  contribution  to  the  observed 
brightness  from  the  Rayleigh  scattering  and  the  contributions  from  the  aerosol  scattering 
and  ground  reflection, 

--  An  uncertainty  in  the  knowledge  of  the  effective  ground  albedo  will  produce  errors 
beyond  our  requirements  in  the  recovered  atmospheric  density  profile. 

Stellar  occultation  and  Resonance  fluorescence  of  O2  Techniques 

•  are  feasible 

-  both  techniques  use  classical  theory  of  photoabsorption  spectroscopy, 

-  in  addition  to  the  atmospheric  density  profile,  both  techniques  can  recover  the  temperature 
profile  since  the  absorption  cross-sections  used  in  both  techniques  are  temperature 
dependent. 

5.2  Accuracy 
5.2.1  Statistical  Error 
Stellar  Occultation  Technique 

•  varies  depending  upon  the  data  analysis  techniques  and  the  instrument  sensitivity, 
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•  varies  from  one  profile  to  another  due  to  the  changes  of  star  brightness, 

•  cannot  be  improved  by  the  integration  time, 

—  the  integration  time  is  fixed  by  the  orbital  velocity  and  the  star  location. 

•  can  be  improved  by  using  barostatic  equation  since  the  temperature  profile  is  also  obtained. 

Resonance  Fluorescence  of  O 2  Technique 

•  also  varies  depending  upon  the  data  analysis  techniques  and  the  instrument  sensitivity, 

•  is  similar  from  one  profile  to  another  due  to  a  steady  emission  brightness, 

•  can  be  improved  by  the  integration  time, 

--  for  example,  the  emission  is  weaker  at  higher  altitudes,  the  accuracies  in  the  recovered 
density  and  temperature  can  be  improved  by  increasing  the  integration  time. 

•  can  also  be  improved  by  using  barostatic  equation. 

5.2.2.  Systematic  Error 
Stellar  Occultation  Technique 

•  arises  due  to  the  uncertainty  in  the  absorption  cross-sections, 

—  4  percent  and  somewhat  greater  for  a  <  5.0  x  1022  cm2  (Yoshino  et  al.,  1987), 

•  arises  due  to  the  uncertainty  in  the  absorption  line  widths, 

—  negligible  considering  the  few  percent  error  in  the  measured  widths  (Yoshino,  1988, 
private  communication) 

•  arises  due  to  the  uncertainty  in  the  line  positions, 

—  negligible  consideming  the  0.0013  nm  accuracy  in  the  measured  line  positions  (Cheung 
et  al.,  1986b) 

•  arises  due  to  the  uncertainty  in  the  temperature  dependency  of  the  cross-sections. 

—  negligible  since  the  energies  of  both  the  upper  and  lower  states  are  accurately  known. 

Resonance  Fluorescence  of  O 2  Technique 

•  arises  due  to  the  uncertainty  in  the  emission  and  absorption  line  strengths, 

—  ~10  percent  and  varies  from  line  to  line  (AFGL  line  compilation  of  McClatchey  et  al., 
1973), 

—  can  be  improved  by  using  more  accurately  measured  line  strengths  (photo-accoustic 
spectroscopy  technique  by  Paul  B.  Hays  at  SPRL,  the  University  of  Michigan,  1988) 
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•  arises  due  to  the  uncertainty  in  the  Franck-Condon  Factor  (ratio  between  the  transition 
probabilities  of  the  O2  Atmospheric  (0-0)  band  and  (0-1 )  band), 

-  may  be  few  tens  of  percent  if  calculated  Franck-Condon  factors  of  Nicholls  (1965)  are 
used, 

--  can  be  improved  if  more  accurate  value  determined  from  this  resonance  fluorescence 
experiment  is  used. 

•  arises  due  to  the  uncertainty  of  the  line  positions 

--  negligible  considering  the  accuracy  in  our  knowledge  of  the  line  positions, 

•  arises  when  the  temperature  dependency  on  the  cross-sections  is  considered, 

--  negligible  since  the  energies  of  both  the  upper  and  lower  states  are  accurately  known. 

5.2 .3.  Pointing  and  Satellite  Position  Error 

•  arises  due  to  the  uncertainly  in  the  tangent  height  determination, 

--  can  be  improved  by  having  accurate  pointing  and  satellite  tracking 

5.3  Complexity 

5.3  1  Operation 

Stellar  Oc  cultation  Technique 

•  requires  accurate  knowledge  of  the  satellite  position  all  the  time  as  planning  aids, 

•  requires  accurate  pointing  in  locating  the  star  as  the  source  of  light, 

•  requires  an  on-board  computer  and  an  efficient  algorithm  to  locate  the  best  available  star  for 
the  next  occultation  sequence  for  autonamous  operation, 

•  involves  frequent  azimuthal  and  zenith  movements  of  the  telescope. 

Resonance  Fluorescence  of  O2  Technique 

•  does  not  require  accurate  knowledge  of  the  satellite  position  all  the  time, 

•  requires  precise  pointing  for  the  accurate  determination  of  tangent  ray  height, 

•  involves  little  azimuthal  movement  of  the  telescope,  only  zenith  scan  is  required. 


5.3.2  Data  Analysis 


Stellar  Occupation  Technique 

•  is  simple  and  efficent 
Resonance  Fluorescence  of  O2  Technique 

•  is  less  efficient  because  it  involves  the  analysis  of  both  the  (0-0)  and  (0-1)  band  emissions. 
5.4  Data  Coverage 

5.4. 1  Spatial  Coverage 
Stellar  Occupation  Technique 

•  covers  an  altitude  region  between  70  and  200  km, 

•  covers  a  latitude  region  determined  by  the  orbital  inclination, 

•  covers  geographical  locations  scattered  around  the  orbital  track, 

•  has  limited  capability'  of  studying  particular  region  of  interest. 

Resonance  Fluorescence  of  O 2  Technique 

•  covers  an  altitude  region  between  60  and  1 50  km, 

•  covers  a  latitude  region  determined  by  the  orbital  inclination, 

•  covers  geographical  locations  systematically  along  the  orbital  track. 

•  has  the  capability  of  studying  particular  region  of  interest. 

5.4.2  Temporal  Coverage 
Stellar  Occupation  Technique 

•  requires  special  consideration  to  the  instrumental  design  in  order  to  have  a  capability  of 
obtaining  a  complete  local  time  coverage  (to  perform  daytime  measurements), 

—  daytime  measurements  require  small  field-ot-view  to  decrease  the  contribution  from  the 
atmospheric  emissions  and  the  scattered  sunlight.  Very  accurate  telescope  pointing  and 
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knowledge  of  satellite  position  thus  are  required. 

•  also  requires  daytime  measurement  capability  in  order  to  be  able  to  obtain  a  complete  seasonal 
coverage  (for  high  inclination  orbits), 

--  there  will  be  periods  of  time  in  a  year  that  the  satellite  is  continuously  exposed  in  the 
twilight  condition  for  high  inclination  orbits. 

•  time  scale  is  not  continous  because  of  the  frequent  telescope  azimuthal  movements, 

•  has  more  dead  time  (time  that  instrument  is  not  taking  measurement). 

Resonance  Fluorescence  of  O2  Technique 

•  has  the  capability  of  obtaining  a  complete  local  time  coverage, 

•  has  the  capability  of  obtaining  a  complete  seasonal  coverage, 

•  time  scale  is  continous, 

•  has  less  dead  time. 

5.5  Scientific  Merits 

Stellar  Occultation  Technique 

•  is  a  dedicated  technique  to  measure  the  density  and  temperature  profiles  of  the  atmosphere. 
Resonance  Fluorescence  of  O 2  Technique 

•  allows  us  to  investigate  the  photochemistry  of  the  O2  Atmospheric  band  emissions, 

•  allows  us  to  examine  the  role  that  the  O2  photodissociation  plays  in  the  Earth's  thermal 
budget  (see  text), 

•  allows  us  to  study  the  effects  of  planetary  and  gravity  waves  in  the  mesospheric  region  not 
only  from  the  measured  perturbations  of  density  and  temperature,  but  also  from  the  airglow 
emissions. 

5.6  Summary 

The  Rayleigh  scattering  technique  is  not  a  viable  technique  to  measure  the  atmospheric  density  in 
the  region  between  60  and  120  km  compared  to  other  two  techniques:  the  stellar  occultation  and  the 
resonance  fluorescence  of  O2.  The  latter  two  are  equally  attractive  and  both  are  worthy  to  be 
considered  and  tested  for  confirmation  of  science  and  accuracy  requirements. 


6.  CONCEPTUAL  INSTRUMENT  DESIGN 


The  conceptual  design  of  an  instrument  for  the  remote  measurement  of  upper  atmospheric  density 
includes  three  basic  components,  i.e.,  a  scanning  telescope,  a  light  dispersion  unit,  and  a  detector 
(Figure  6.1).  Two  instrument  types  are  the  leading  candidates  for  the  trade-off  study  during  this 
conceptual  design  phase;  the  Image  Plane  Detector  (IPD)  Spectrophotometer  and  the  grating  spec¬ 
trometer  (GS). 

6.1  Image  Plane  Detector  Spectrophotometer  (IPDS) 

The  IPDS  is  an  evolutionary  combination  of  a  variety  of  photometer  using  the  Image  Plane 
Detector  developed  earlier  for  satellite  instruments  (Hays,  et  al.,1981).  It  is  very  similar  to  a  tilting 
filter  photometer.  It  uses  a  multichannel  Image  Plane  Detector  as  the  spectral  selector  to  detect  the 
emission  spectrum  dispersed  by  an  interference  filter.  An  extensive  description  of  this  low 
resolution  instrument  is  included  in  the  appendix,  and  only  a  brief  description  is  given  here. 


FIGURE  6.1 
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MAGE  PLANE  DETECTOR  SPtsCTRGPHOTCMSTEH 
Figure  6.2 

Figure  6.2  presents  a  schematic  diagram  of  the  instrument.  Energy  collection  is  provided  by  a 
telescope  coupled  to  a  scan  mirror  (or  a  scanned  telescope)  with  a  Fibre  optic  line-to-circle 
convener.  The  line-to-circle  converter  transforms  the  rectangular  altitude  elements;  narrow  vertical 
resolution  and  wide  horizontal  resolution,  to  a  circular  output.  The  circular  output  of  the  energy 
collection  subsystem  then  feeds  the  optical  energy  to  a  wavelength  determination  subsystem  which 
consists  of  a  dielectric  interference  filter.  The  interference  filter  has  a  property  that  the  wavelength 
of  maximum  transmission  varies  with  the  incident  angle  of  incoming  light.  It  therefore  disperses 
the  light  coming  from  all  angles  and  generates  an  interference  pattern  which  is  imaged  by  the 
objective  lens  onto  the  Image  Plane  Detector  (1PD). 

The  function  of  the  IPD  here  is  twofold.  First,  it  converts  light  to  electric  signals.  Light  energy  is 
converted  to  photoelectrons  at  the  photo-cathode,  and  proximity-focussed  to  the  input  of  the  micro- 
channel  plate  imaging  multipliers.  The  imaging  electric  multipliers  amplify  the  signals  by  ~106  and 
collect  then  in  different  concentric,  constant  area  metallic  anodes.  The  second  function  of  an  IPD  is 
to  collect  the  light  dispersed  by  the  filter  so  that  each  anode  samples  a  different  pan  of  the  spectra 
with  equal  wavelength  interval.  The  signals  from  the  anodes  are  then  fed  to  the  amplifiers  and  sent 
out  to  the  data  collection  system  for  data  processing. 
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6.2  Ebert-Fastie  Grating  Spectrometer  (EFGS) 


Two  experiments  meet  the  criteria  for  the  measurement  of  neutral  density  in  the  altitude  region  60- 
120  km,  i.e.,  stellar  occultation  and  O2  resonance.  In  the  previous  section  we  discussed  the 
conceptual  design  of  IPDS  as  it  could  be  used  for  the  remote  measurement  of  atmospheric  density 
by  looking  at  spectral  features  of  the  Atmospheric  bands  of  molecular  oxygen,  i.e.,  O2  resonance. 

In  this  section  we  discuss  the  conceptual  design  of  the  Ebert-Fastie  Grating  Spectrometer  as  it 
could  be  used  for  the  measurement  of  O2  resonance  or  stellar  occultation.  Both  types  of  EFGS 
configurations  consist  of  l/8th  meter  Ebert-Fastie  monochromators  with  off-axis  telescopes 
containing  spherical  mirrors.  Both  contain  plane  gratings  to  disperse  the  image  over  the  desired 
spectral  range  and  both  detectors  would  be  Charge-Coupled-Devices  (CCD)  (Figure  6.3). 
Although  the  monochromator  is  expected  to  be  similar  for  either  measurement  scheme  the  design  of 
the  telescope,  the  plane  grating,  and  the  CCD  are  unique  to  the  individual  measurement  technique 
as  are  the  operational  requirements. 

1/8  METER  SPECTROMETER 
SYSTEM 


G.  1/8  METER  SPECTROMETER 

H.  CCD  DETECTOR 


Figure  6 . 3 
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6.2. 1  The  O2  Resonance  Measurement 

6.2. 1.1  Field-of-View  Determination 

It  is  desirable  in  a  limb  measurement  to  use  a  low  altitude  satellite  to  obtain  a  large  field  of  view  for 
a  given  spatial  resolution.  For  the  purpose  of  estimating  the  requirements  on  the  instrument  that 
will  be  used  for  these  measurements  we  can  choose  a  satellite  with  a  circular  orbit  of  550  km.  The 
orbital  inclination  angle  is  arbitrary  at  this  point  and  depends  upon  the  decision  to  concentrate  on 
the  coverage  of  the  earths  surface  (high  inclination)  or  concentrate  on  coverage  of  the  local  times 
(low  inclination).  For  a  satellite  positioi  ed  550  km  above  the  earth  the  line  of  sight  distances  to  the 
tangent  points  at  altitudes  of  60  km  and  1 20  km  are  approximately  2560  km  and  2400  km 
respectively.  If  a  height  resolution  is  selected  to  be  2  km,  the  field  of  view,  in  the  vertical,  will  be 

fy  1cm  1 

8  -  “ilSotetJ  -  °044938” 

=  2.697  arc  min.  (6.1) 

6.2. 1.2  Telescope  Selection  and  Input  Slit  Width  Determination 

The  initial  selection  for  me  telescope  has  a  focal  length  of  250  mm  with  an  F  number  of  f/5.  If  we 
match  the  telescope  req  luements  with  the  required  field  of  view  (F.O.V.)  for  the  2  km  resolution 
we  determine  the  entrance  slit  width  to  be  w=0. 196  mm,  i.e., 

w  =  (250  mm)  tar  (0.044938°)  (6.2) 

This  value  relates  to  a  tangent  altitude  of  60  km.  An  equivalent  calculation  shows  that  the  altitude 
resolution  at  120  km,  ur  ng  this  slit  dimension,  would  be  1.88  km.  We  should  note,  also,  that  we 
would  need  to  install  t.  simple  light  filter  to  prevent  interference  from  low  order  harmonic 
frequencies. 

The  numerical  aperture  for  the  telescope,  and  spectrometer,  is  calculated  to  be 

f/5spec  =  f/5tei  =  1/2  N.A.  =  l/2sin  0,  (6.3) 

2sin  0  =  1/5  (6.4) 
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0  =  sin1  (0.1)  =  5.739  degrees  (half  angle) 


(6.5) 


» 
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Thus  the  total  field  of  view  for  the  telescope  would  be  1 1.478  degrees. 

6. 2. 1.3  Grating  Selection 

In  order  to  maximize  the  signal  level  for  our  measurement  we  choose  an  Ebert-Fastie  spectrometer 
with  a  focal  length  of  125  mm  (1/8  meters).  For  a  spectrometer  of  f/5  the  grating  size  is  given  by 

Dg  cos  0  =  125/5  mm  =  25  mm.  (6.6) 

where  0  is  the  angle  of  rotation  as  determined  by  the  spectrometer  axis  and  the  normal  vector  to  the 
plane  of  the  grating.  This  angle  will  be  determined  in  a  later  section. 


EBERT-FASTIE  GRATING  SPECTROMETER 


Focal  Plane  Eben  Mirror 
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The  grating  equation  is  given  as 


d[sin(0i)  +  sin(0e)]  =  mX,  (6.7) 

where 

0,  =  incidence  angle  of  the  light  ray, 

0e  =  exit  angle  of  the  light  ray, 
d  =  separation  of  the  ruled  lines, 
m  =  interference  order  number, 

X  =  wavelength  of  the  light  ray. 

The  angles  0*  and  0e  are  measured  with  respect  to  the  coordinates  of  the  plane  grating  where  the  x- 
y  plane  coincides  with  the  grating  surface  and  and  the  z  axis  is  normal  to  the  surface.  These  angles 
can  be  measured  with  respect  to  the  spectrometer  coordinate  system,  as  shown  in  Fig.  6.4  where 

51  =  inlet  slit, 

52  =  exit  slit  (position  of  the  detector), 

0  =  the  angle  defined  by  the  z-axis  of  the  plane  gTating  and  the  axis  of  the  spectrometer, 

a  =  the  angle  determined  by  the  incident  and  exit  light  rays  with  respect  to  the  spectrometer 
axis.  Note:  this  angle  is  fixed  and  is  determined  only  by  the  spectrometer  geometry. 

The  incident  and  exit  angles  can  be  redefined  for  the  spectrometer  coordinate  system  as 

0(  =  0  -  a  (6.8) 

0e  =  0  +  ot 

and  grating  equation  becomes 
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d[sin(9  -  a)  +  sin(0  +  a)]  =  mX, 


(6.9) 


or 

2d  sin  0  cos  a  =  mX.  (6.10) 

This  equation  is  specific  for  the  Ebert-Fastie  spectrometer  configuration.  The  angle  a  is  a  function 
of  the  geometry,  but  for  the  purpose  of  preliminary  design,  may  be  assumed  to  have  a  value  on  the 
order  of  8- 10  degrees.  The  quantities  d,  m,  and  0  all  depend  on  the  wavelength,  X,  and  in  general 
are  a  compromise  depending  on  the  specific  instrument  requirements. 

Grating  requirements 

The  wavelength  region  for  this  measurement  lies  in  the  near  infrared  (N1R),  7650A  and  8650A.  If 
we  select  7650A  for  X,  and  for  a  starting  point  choose  d  =  1/1800  lines/mm,  e.g.,  these  are  the 
values  that  the  NUV  spectrometer  uses  in  the  RAIDS  (McCoy,  et  al„  1986)  experiment,  then  we 
can  determine  the  values  of  0  as  a  function  of  the  interference  order  by  using  Eq.  (6. 10). 

For  example,  at  m  =  1,  and  X  =  7650  A 

0  =  sin'1  [ Y^cofa  ]  <6-ll> 

=  44.193°. 

If  we  repeated  the  last  calculation  but  used  a  X  equal  to  8650A  we  could  obtain  9  equal  to  52.018°. 
Physical  size 

Using  the  9  of  44.193°  required  for  a  wavelength  at  7650A,  we  can  calculate  the  dimension  of  the 
grating,  Dg,  for  a  125  mm  spectrometer  at  f/5  from  Eq.  (6.6), 

Dg  =  30.684  mm. 

Repeating  this  calculation  for  X  equal  to  8650A,  where  9  is  52.018°,  the  value  of  Dg  becomes 
40.623  mm. 

If  we  are  to  maximize  the  signal  at  8650A  we  probably  should  use  the  larger  physical  size  for  the 
grating. 
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Grating  Dispersion 


The  dispersion,  Pe,  for  a  grating  spectrometer  is  calculated  by 

Pe  =  (mfgfr  )/d  cos(6e )  (6.12) 

where 


feff.  =  125  mm, 

9e=  53.193°  (at  A  equal  to  7650  A), 
d  =  1/1800  mm, 
m  =  1. 

The  dispersion  at  the  exit  position  in  the  focal  plane  is  then  calculated  to  be 
Pe  =  0.037555  mm/A, 
or  about  26  6276 A/mm. 

Thus  the  resolution  fora  slit  width  of  0.196  mm  would  be  equal  to  5.219A. 

6. 2. 1.4  Detector  Selection:  A  CCD  Device 

First,  we  need  to  assume  something  about  the  CCD.  For  example,  we  need  to  know  something 
about  its  array  size,  its  dark  current,  its  read  noise,  its  spectral  range,  its  charge  collection 
efficiency,  the  pixels  sizes,  and  probably  some  other  things  but  this  should  get  things  started.  Lets 
assume  that  we  have  a  Tektronix  CCD  (Janesick.  et  al.,  1987).  Its  array  size  is  512  x  512  and  the 
indvidual  pixel  size  is  27  um.  The  physical  array  size  is  13.8  x  13.8  mm2.  A  CCD  of  this  physical 
size  gives  a  length  dimension  of  the  inlet  slit  no  greater  than  13.8  mm  long.  It  also  determines  the 
required  dispersion.  For  example,  if  we  w  ish  to  disperse  the  emission  band  over  the  entire  face  of 
the  CCD  we  need  a  large  dispersion.  Using  the  dispersion  of  the  exit  ray  at  the  focal  plane  of 
0.037555  mm/A,  or  inversely  26.6276A/mm,  one  would  obtain  a  spread  of  367. 5A  over  the  face 
of  the  CCD. 


Carrying  out  the  calculation  of  the  dispersion  for  X  =  8650  A,  we  get 
Pe  =  0.046436  mm/ A, 

or  inversely,  21.53488A/mm.  A  wavelength  spread  of  291  ,lk  over  the  13.8  mm  width  of  the 
CCD  is  then  obtained. 

The  two  O2  Atmospheric  bands  that  we  are  considering  both  cover  a  spectral  region  about  200A.  It 
is  therefore  reasonable  at  this  point  to  consider  the  initial  design  as  adequate. 


6.2.2  The  Occultation  Measurement 

6.2.2. 1  Telescope  Selection 

First  we  start  with  an  f/5  125  mm  Ebert-Fastie  spectrometer.  The  telescope  should  be  matched  to 
the  spectrometer,  i.e.,  it  should  have  the  same  F  number.  Using  an  objective  lens  diameter  of  40 
mm  which  we  run  a  statistical  error  analysis  on  (See  Section  2),  we  obtain,  for  an  f/5  telescope,  a 
corresponding  focal  length  of  200  mm  and  numerical  aperture  (N.A.)  of  5.739  degrees.  Since  the 
total  Field  of  view  of  11.478  degrees  (  2  x  N.A.)  would  be  too  large  for  the  occulation 
measurement,  a  field  stop  is  required  in  order  to  to  exclude  the  stars  other  than  the  one  selected. 
The  field  stop  would  be  the  effective  input  slit  for  the  spectrometer. 

6.2. 2. 2  Grating  Selection 

In  order  to  maximize  the  signal  level  for  this  measurement  we  need  to  match  the  grating  size  with 
the  spectrometer,  i.e.,  for  an  f/5  spectrometer  with  a  focal  length  of  125  mm, 

DgCosQ  =  (125/5)  mm  =  25  mm  (6.13) 

where  9  is  the  angle  of  rotation  as  determined  by  the  spectrometer  axis  and  the  normal  vector  to  the 
plane  of  the  grating  surface,  and  will  be  determined  later. 

The  grating  equation  is  given  as 

dfsin  0j+ sin  0el  =  mX,  (6.14) 
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where 


6j  =  0  -  a 
0e  =  0  +  a 

and  equation  (6.14)  can  be  reduced  to  give 

2d  sin  0  cos  a  =  mX.  (6. 1 5) 

Note  that  a  is  a  constant  angle  as  determined  by  the  geometry  of  the  spectrometer.  Its  value  is  on 
the  order  of  8°-10°.  We  choose  a  value  of  9°  for  our  calculations. 


Grating  Requirements 

The  wavelength  region  for  this  measurement  lies  in  the  Schumann-Runge  Bands,  1750-2000  A.  If 
we  select  1750  A  for  X,  and  for  a  starting  point  choose  d  =  1/1200  lines/mm  we  can  determine 
values  of  0  as  a  function  of  the  interference  order  by  using  equation  (2). 

For  example,  at  mp=  1,2,..,  and  X  =  1750  A  then 


0  =  sin  1 


m  -  1750 
2d  cos  a 


(6.16) 


and  we  get  the  values  in  Table  6. 1 . 


Table  6.1 


m 

e 

1 

6.10° 

2 

12.28° 

3 

18.60° 

4 

25.17° 

5 

32.11° 

6 

39.63° 

7 

48.09° 
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Physical  Size 


Using  the  6  of  18.60°  at  m=3  required  for  a  wavelength  at  1750A,  we  can  calculate  the  dimension 
of  the  grating,  Dg,  for  a  1 25  mm  spectrometer  at  fJ5  from  Eq.  (6.6), 

Dg  =  26.4  mm. 

Grating  Dispersion 

The  dispersion,  Pe,  for  a  grating  spectrometer  can  be  calculated  by 

Pe  =  (mfeff.)/d  cos(8e)  (6.17) 

where 

feff.  =  125  mm, 
a  =  9°  (  assumed  ), 
m  =  3, 

ee=  18.60°  +  9°  =  27.60°, 
d  =(1/1200)  mm. 

The  dispersion  at  the  exit  position  in  the  focal  plane  of  the  spectrometer  is  then 
Pc  =  0.05078  mm/A, 
or  19.693  A/mm. 

6. 2.2. 3  Detector  Selection:  A  CCD  Device 

The  CCD  assumed  for  our  instrument  has  an  array  size  of  512  x  512  pixels  square.  The  physical 
array  size  is  13.8  x  13.8  mm2.  For  a  dispersion  0.05078  mm/A  we  calculate  that  we  would  be  able 
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to  spread  271  A  over  the  surface  of  the  CCD,  i.e.,  we  would  be  able  to  record  on  the  CCD  the 
entire  wavelength  spectrum  for  each  selected  integration  period.  The  integration  period  would 
corres-pond  to  the  period  of  time  for  the  tangent  height  ray  to  change  by  2  km  during  the 
occultauon  sequence  as  described  earlier. 


6.3  Summary  of  Instrument  Parameters  for  EFGS 


1 

Ch  Resonance 

Stellar  Occultation 

7550,  8750 

1750-2000 

F-number 

f/5 

f/5 

Objective  Diameter  (mm) 

50.0 

40.0 

Field  of  View  (deg2) 

0.045  x  6.7 

3.0  x  3.0 

Input  Slit  Size  (mm2) 

0.196  x  29.1 

10.5  x  10.5 

Grating  Area  (mm2) 

40  x  40 

27  x  27 

Grating  Density  (lines/mm) 

1800 

1200 

Interference  Order 

1 

3 

Exit  Aperture  (mm2) 

13.8  x  13.8 

13.8  x  13.8 

Detector  Array  Size 

512  x  512 

Spectral  Coverage  (A) 

368,  297* 

250 

Detector  Pixel  Size 

27  x  27  um2 

27  x  27  um2 

the  spectral  coverage  is  given  for  the  lines  at  7650  A  and  8650  A  respectively. 


7.  OPERATIONAL  REQUIREMENTS 


Although  some  of  the  operational  requirements  of  the  O2  resonance  and  stellar  occultatic 
measurements  are  similar,  fundamentally  the  requirements  are  quite  diverse.  The  optimum  ort 
plane  for  either  measurement  is  expected  to  have  an  inclination  on  the  order  of  60°.  This  inclii 
ation  angle  would  provide  global  coverage  of  the  earth's  surface  while,  at  the  same  time,  providir 
reasonable  local  time  coverage.  It  is  also  desirable  for  the  satellite  to  have  a  relatively  low  ear 
orbit  for  either  measurement  technique  but  for  different  reasons. 

7.1.  02  Resonance 

There  are  two  reasons  why  it  is  desirable  to  have  a  low  earth  orbit  for  the  O2  resonanc 
measurement.  The  2  km  height  resolution  imposed  for  the  density  measurement  places  a  restrictic 
on  the  spectrometer  input  slit  width  dimension  that  is  inversely  proportional  to  the  distance  betwee 
the  position  of  the  satellite  and  the  tangent  height  altitude.  Any  reduction  in  slit  width  is  inverse! 
proportional  to  the  amount  of  light  flux  that  can  be  passed  through  the  instrument  which  effective! 
reduces  the  signal  to  noise.  But,  perhaps,  a  more  important  consideration  relates  to  the  pointin 
requirements  for  resonance  measurement.  We  not  only  need  to  resolve  2  km  altitude  increment 
discussed  above,  over  the  altitude  range  60  -  120  km  but  we  need  to  know  precisely  where  the: 
height  increments  lie  such  that  we  can  calculate  the  integrated  contents  of  each  height  packaj 
within  5%.  If  we  assume  an  isothermal  atmosphere  then  the  density  can  be  expressed  as 

\h 

p  =  p,  e  11  (7.1) 


where 

p0  =  density  at  h„, 
p  =  density  at  h0  +  Ah, 

H  =  density  scale  height. 

The  uncertainty  in  the  competed  density  due  only  to  the  uncertainties  in  these  parameters  can 
expressed  as 


^-dnTjiE. 

,riu  :  i  j  .  o.\h 


dDh! 


6p  = 


(7.2 


where 


\j/p  =  28.2  seconds'1. 


Thus  for  a  remote  measurement  of  O2  resonance  the  pointing  requirements  are  rather  stringent. 
The  pointing  accuracy  would  have  to  be  on  the  order  of  20  arc  seconds. 

While  the  pointing  accuracy  is  rather  severe  for  this  measurement  technique  accurate  knowledge  of 
the  satellite  position  is  not  critical.  The  uncertainty  in  this  parameter  is  translated  into  an 
uncertainty  in  the  measurement  track.  The  measurement  bins,  i.e.,  the  horizontal  dimensions,  are 
not  expected  to  be  impacted  by  normal  tracking  uncertainties. 

Both  of  the  instrument  types  that  could  be  employed  for  the  O2  resonance  measurement,  the  Image 
Plane  Detector  Spectrophotometer  (IPDS)  and  the  Ebert-Fastie  Grating  Spectrometer  (EFGS),  have 
the  same  satellite  operational  requirements.  In  addition  to  the  pointing  requirements  described 
above,  the  instrument  telescope  must  be  capable  of  360  degree  movement  in  azimuth  and  approxi¬ 
mately  10  degrees  in  elevation  with  respect  to  a  stabilized  platform  that  is  fixed  horizontal  to  its 
sub-satellite  position.  The  instrument  telescope  should,  at  all  times,  be  pointed  at  -90°  away  front 
the  sun  vector  to  avoid  the  scattered  sunlight.  The  measurement  can  be  obtained  in  daytime  and/or 
nighttime  conditions. 

7.2  Stellar  Occupation 

As  we  have  discussed  in  the  previous  section  on  the  O;  resonance  technique,  the  atmospheric 
density  was  resolved  in  2  km  layers  by  imaging  the  contents  of  a  spectrometer  input  slit,  of 
appropriate  width,  onto  a  CCD  while  the  telescope,  or  telescope  mirror,  was  sequencially  scanned 
over  the  60-120  km  altitude  region.  We  also  found  that  the  accuracy  in  the  measured  density 
profile  w  as  tied  directly  to  our  knowledge  of  the  pointing  direction  and  that  if  we  want  to  measure 
w  ith  5%  accuracy  we  need  good  pointing  accuracy,  on  the  order  of  20  arc  seconds. 

In  contrast,  the  stellar  occultation  measurement  resolves  the  layered  density  information  by 
selecting  appropriate  integration  periods  as  the  occulted  star  descends  on  the  earth's  horizon  aft  of 
the  satellite  motion.  It  is  desirable  to  have  a  low  orbital  altitude  to  increase  the  occultation  period 
and  corresponding  integration  time  per  altitude  increment. 

Also  in  section  7.1  we  found  the  necessity  of  resolving  the  tangent  height  to  approximately  0.35( 
km  in  order  to  obtain  a  5%  measure  of  atmospheric  density.  For  the  occultation  measurement  ar 
error  in  satellite  position  induces  a  tangent  height  measurement  error.  If  we  refer  to  Figure  7. 1  wi 
see  that 


Figure  7.1 


cos  (O)  =  (re  +  h)  /  (re  +  hs) 
h  =  (re  +  hs)  cos  (Q)  +  re 


where 

W  =  satellite  nadir  angle  =  21.69°  (at  assumed)  condition), 
re  =  radius  of  earth  (assumed  to  be  637 1  km), 
hs  =  satellite  altitude  (assumed  to  be  550  km), 
h  =  tangent  height  altitude. 

Ws  =  satellite  radian  frequency 

=  1.0957j6t-3  rad./sec.  (at  assumed  altitude) 

Note:  Off  orbital  track  error  considered  negligible  for  this  calculation. 


(7.9) 

(7.10) 


The  uncertainty  in  the  tangent  height  altitude  as  a  function  of  the  variable  parameters  can  be 
expressed  as 


js-fL  f  (|Lfrdwf 
^  J  {  I  i  j  ■ 


(7.11) 


where 


~  cos  £i. 
oh,. 


=  -  (re  +  hs)  sin  ft. 


(7.13) 


8  ft  =  Ws5t. 


(7.14) 


From  equation  (7.6)  we  see  that,  in  order  to  resolve  0.350  km  in  the  tangent  height  altitude  the 
uncertainty  in  the  satellite  altitude  would  have  to  be 


5h,  <  -  0  377  km 

cos  ft  cos  21.69 


(7.15) 


and  the  uncertainty  in  the  nadir  angle  would  have  to  be 


(r.  +  hs)  sin  0) 


(7.16) 


•e  1 .368  x  10  radian. 


and  we  can  translate  the  uncertainty  in  angle  to  an  uncertainty  in  time,  through  Eq.  (7.7), 


5t  <  1  -3-6g-x  ML—  =  0.125  seconds 


(7.17) 


1.095755  x  10 


These  calculations  were  made  for  a  specific  satellite  altitude  and  tangent  height  measurement  lr 
order  to  establish  the  order  of  magnitude  for  the  parameter  uncertainties. 
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Abstract 


A  new  variety  of  low  resolution  spectrometer  is  described  in  this  study.  This  device,  an  Image 
Plane  Detector  Spectrophotometer,  has  high  sensitivity  and  modest  resolution  sufficient  to 
determine  the  rotational  temperature  and  brightness  of  molecular  band  emissions.  It  uses 
interference  filter  as  dispersive  element  and  a  multichannel  Image  Plane  Detector  (IPD)  as  photon 
collecting  device.  The  data  analysis  technique  used  to  recover  the  temperature  of  the  emitter  ?nd  the 
emission  brightness  is  presented.  The  Atmospheric'  band  of  molecular  oxygen  is  used  to  illustrate 
the  application  of  the  device. 
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I.  Introduction 


The  temperature  of  emitting  atoms  or  molecules  is  reflected  in  their  spectral  line  shape  or  the 
emission  band  structure.1  The  width  of  an  individual  atomic  or  molecular  line  can  be  used  to  obtain 
the  doppler  temperature.2,3  This  is,  however,  difficult  for  the  emissions  originating  in  the 
mesosphere  where  the  temperature  is  low  and  the  line  width  is  narrow.  Under  these  circumstances, 
a  high  resolution  instrument  is  required.  An  alternative  remote  sensing  technique  to  monitor  the 
temperature  is  to  obtain  the  rotational  temperature  from  the  molecular  band  structure  by  using  a 
scanning  spectrometer.4,5  A  finite  period  of  time  is  usually  required  to  scan  a  whole  emission  band 
system,  consequently  lowering  the  temporal  resolution.  For  a  measurement  made  by  instruments 
on  board  rocket  or  satellite,  this  would  result  in  a  lower  spatial  resolution.  Our  intent  here  is  to 
design  a  spectrophotometer  with  a  relatively  high  throughput  and  with  spectral  resolution  sufficient 
to  resolve  the  molecular  band  structure  w  ithout  the  necessity  of  scanning. 

The  tilting-filter  photometer  has  been  used  frequently  in  the  past  to  examine  the  emissions  in  the 
airglow  and  aurora.6'8  It  has  a  spatially  fixed  aperture,  and  the  wavelength  scan  is  carried  out  by 
mechanically  tilting  the  filter  with  respect  to  the  optical  axis  of  the  instrument.  Although  it  has 
relatively  high  throughput,  its  instrument  broadening  varies  with  tilt  angles  of  the  filter,  and 
consequently,  the  wavelength.  In  order  to  correct  this  problem,  a  multi-spiral  sector 
spectrophotometer  was  developed.4  In  this  instrument  a  field  stop,  made  of  a  sector  shaped 
masking  disk  and  a  disk  with  a  spiral  slit,  was  put  in  front  of  the  detector  so  that  the  signals  of  the 
emission  can  be  selectively  collected  at  a  given  wavelength  range.  The  spectral  scan  is  accomplished 
by  rotating  the  spiral  slit,  and  is  capable  of  obtaining  a  nearly  constant  peak  transmittance  and  band 
pass  for  large  tilt  angles.  This  instrument  has  been  successfully  used  in  monitoring  auroral  0+  (~P  - 
2D)  emissions  at  7320A  by  Swenson  ei  al . 4  Both  the  tilting-filter  photometer  and  the  spiral 
spectrophotometer  scan  the  selected  spectrum  by  either  tilting  the  interference  filter  to  shift  the  center 
wavelength  or  moving  a  shaped  aperture  across  the  interference  pattern  created  by  a  dielectric 
interference  filter,  while  the  detector  (photomultiplier)  is  in  a  fixed  position. 

The  instrument  described  here  is  very  similar  in  principle  to  the  tilting-filter  photometer  and  the 
spiral  spectrophotometer.  It  uses  a  multichannel  Image  Plane  Detector  (1PD)  as  the  spectral  selector 
to  detect  the  emission  spectrum  dispersed  by  an  interference  filter.  The  Image  Plane  Detector  was 
first  used  on  the  Fabry-Perot  Interferometer  flown  on  the  Dynamic  Explorer  Satellite  10  The 
geometrical  parameters  and  characteristics  of  the  IPD  have  been  described  in  details  by  Killeen  e 
al.n  In  brief,  the  IPD  is  an  electro  optical  device  which  converts  the  intensity  distribution  of  ligh 
falling  on  the  image  plane  into  a  set  of  discrete  electron  pulse.  It  collects  the  light  falling  on  severa 
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equal  area  annular  anodes,  arranged  in  the  shape  of  a  'Bullseye'  pattern.  This  pattern,  as  we  shall 
see,  is  such  that  every  channel  samples  a  range  of  the  spectrum  with  equal  wavelength  interval.  The 
usage  of  the  IPD  make  this  photometer  not  only  have  almost  identical  transmission  function  for  each 
channel,  but  also  be  able  to  obtain  spectral  information  with  no  mechanical  scanning. 

This  paper  will  describe  a  low  resolution  spectrophotometer  using  Image  Plane  Detector  (IPD) 
as  the  spectral  selector.  The  data  analysis  technique  to  recover  the  brightness  of  the  emission  source 
and  the  rotational  temperature  of  the  emission  band  will  be  discussed  as  will  the  associated  errors. 
The  02  (0-1)  atmospheric  band  emission  will  be  used  as  an  example  for  demonstration. 

FI.  Description  of  the  Image  Plane  Detector  Spectrophotometer 

The  instrument  discussed  here  is  intended  to  measure  the  molecular  emission  band  structure 
with  relatively  low  spectral  resolution.  It  consists  of  three  major  elements:  interference  filter, 
objective  Uns  and  Image  Plane  Detector.  The  schematic  diagram  of  the  instrument  is  shown  in 
figure  1.  The  light  coming  from  a  radiative  source  is  collected  and  passes  through  the  interference 
filter  with  a  wide  angular  range.  The  interference  filter  here  acts  as  a  spectrum  separating  elements, 
generating  an  interference  pattern  which  is  imaged  by  the  objective  lens  onto  an  Image  Plane 
Detector  (IPD).  In  principle,  the  interference  filter  is  a  Fabry-Perot  etalon  with  very  narrow  spacer. 
For  a  monochromatic  incident  light,  a  concentric  Fabry-Perot  ring  pattern,  characterized  by  a 
gradually  narrowing  fringe  width  with  increasing  radial  distance  from  the  optical  a  is,  can  be 
formed.  A  single  fringe  here  corresponds  to  the  maximum  light  transmission  at  the  special 
wavelength  of  the  incoming  light.  The  Image  Plane  Detector  uses  a  concentric-ring  multiple  anode 
and  incorporates  a  three-stage  microchannel  plate  electron  multiplier.  The  unique  anode  patten 
shape  of  the  IPD  was  designed  to  match  the  interference  patten  formed  by  the  filter  or  etalon  so  that 
the  equal  area  anodes  monitor  the  equal  interv  al  of  the  wavelength  spectrum. 

The  relationship  between  incident  angle  of  light  0  to  the  filter  surface  normal  and  the  wavelength 
of  maximum  transmission  X0  is  giv  en  bv  h 


where  X0  is  the  wavelength  of  maximum  transmittance  for  0  =  0°.  and  q*  is  an  effective  refractive 
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index  of  the  filter  which  lies  between  the  two  refractive  indices  of  the  dielectric  materials.  The 
center  wavelength  Xq,  the  width  of  the  filter  and  the  focal  length  of  the  objective  lens  have  to  be 
selected  based  on  the  spectral  coverage  of  the  emission  of  interested.  It  can  be  seen  from  equation 
(1)  that  as  the  incident  angle  of  light  increases,  the  shifts  toward  shorter  wavelength.  Therefore 
in  order  to  use  the  filter  as  the  spectrum  dispersive  element,  the  center  wavelength  of  the  filter  X0 
has  to  be  selected  at  the  longest  possible  wavelength  of  any  emission  band  to  be  monitored.  The 
outmost  anode  of  the  IPD  measures  the  light  coming  from  the  maximum  incident  angle  9max, 
consequendy  determining  the  instrument  field  of  view  and  the  spectral  coverage.  The  width  of  the 
filter  and  the  focal  length  of  the  objective  lens  have  to  be  chosen  so  that  the  instrument  transmission 
function  not  only  be  able  to  cover  the  whole  emitting  band,  but  also  be  able  to  resolve  the  features 
observed  by  two  adjacent  anodes. 

In  order  to  illustrate  this  issue  we  consider  a  theoretical  double  cavity  interference  filter  shown 
in  figure  2,  built  by  following  Lissberger's  theory.12  The  two  refractive  indices  for  the  dielectric 
materials  are  2.30  and  1.35  respectively.  The  center  wavelength,  at  which  the  transmission  is 
maximum  for  normal  incident  light,  is  selected  to  be  8730  A,  corresponding  to  the  longest 
wavelength  of  02  (0-1)  atmospheric  band,  and  the  full  width  at  half  maximum  transmittance 
(FWHH)  is  about  12  A.  The  effective  refractive  index  |i*  for  this  filter  is  about  1.49.  From 
equation  (1),  the  center  wavelength  shift.  Xe  -  Xq,  is  proportional  to  the  square  of  the  incident  angle 
Q2.  Therefore  in  order  for  the  instrument  to  have  a  spectral  coverage  of  about  150  A  (the  spectrum 
range  of  O,  (0-1)  atmospheric  band),  the  maximum  0  should  be  in  the  neighborhood  of  15°. 
Figure  3  shows  the  filter  transmission  functions  Tk  for  twelve  different  incident  angles  0k,  where 

0k  =  (  k  -  1)1/2  0c.  and  k  =  1,  2 .  12.  The  angle  0,  corresponds  to  the  normal  incident  angle 

and  equals  to  zero,  and  angle  02  =  0c  is  selected  so  that  by  using  the  filter  shown  in  figure  2.  the 
transmission  function  at  two  adjacent  angles  0k  and  0k+)  can  be  clearly  resolved.  Here.  0c  is  equal 
to  4.4°.  The  anode  geometry  of  the  detector  allows  us  to  collect  light  signals  at  different 
wavelength  dispersed  by  the  filter,  and  subsequently  defines  the  focal  length  of  the  objective  lens 
required  for  the  instrument.  Note  that  a  1 2  channel  IPD  similar  to  the  one  used  in  the  Fabry-Perot 
interferometer  on  board  the  Dynamic  Hxplorer  Satellite11  is  used  here.  As  the  geometric  parameter 
or  the  number  of  anodes  of  the  IPD  changes.  9c  w  ill  be  different.  The  focal  length  of  the  objective 
lens  in  this  case  is  determined  to  be  2.27  cm. 

Due  to  the  finite  width  of  the  detector  ring,  the  instrument  transmission  function  Tk  at  a  given 
anode  k  is  broadened  and  can  be  evaluated  as 
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(2) 


with 


2k  it 

—  J  d$  J  T(X,0)  Dk(0)  sin  0  d0, 

^  o  o 


Dk(0)  =  1 


0k<  0  <  0k+  A0k 


=  0 


else. 


where  £2k  is  the  field  of  view  of  the  k-th  anode,  T(X,0)  is  the  filter  transmission  function,  and  0k, 
9k+A9k  are  the  angles  defined  by  k-th  anode  of  the  detector  assuming  that  Dk  is  an  ideal  square 
function.  Here  Tk  is  normalized  by  so  that  the  peak  transmission  is  unity. 

For  small  angles  of  0,  equation  (2)  can  be  approximated  as 


Tl(X) 


T(X,8)  Dk(6)  0  d9. 


(3) 


Writing  equation  (3)  in  terms  of  A.  from  its  relationship  with  0,  one  obtains 


?k(M  =  -2±ll  f  T(X,A )  Dk(X  J  dX  (4) 

*  oiJ  8  a  0 

H  A0  o 

The  normalized  square  function  Dk(Ag)  limits  the  integration  between  two  center  wavelengths  which 
correspond  to  the  inner  and  outer  edges  of  the  k-th  detector  ring.  Mathematically,  Tk  can  be 
considered  as  a  convolution  of  the  filter  transmission  function  T(X.0)  and  a  detector  broadening 
square  function  Dk', 


where 


and 


T  =  T  *  D, 


D  '  =  1/AA 
k  ak 

=  0 


U  -a  I 

9 


<  AX 


0k 


else. 


(5) 
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Figure  4  presents  the  simulated  instrument  transmission  functions  Tk  multiplying  by  ftj,  for  the 
12  channels  of  the  IPD.  The  result  is  calculated  based  upon  the  convolution  of  the  double  cavity 
interference  Filter  transmission  function  (  Figure  3)  and  the  transmission  function,  represented  by  a 
square  function,  of  the  12  channel  IPD.  We  note  the  constant  spectral  resolution  for  this  device, 
compared  to  devices  like  tilting  filter  photometer  which  exhibits  large  changes  in  resolution  with 
wavelength  shift. 

If  the  radiative  source  of  interest  is  a  molecular  emission  band,  one  can  obtain  the  output  in 
counts  for  the  k-th  channel  as 

in6  m6  r 

Nk  -  —  R,  A„Q  *  I  B(X)  •  yy  +  —  Rc  A„Q  At  j  T„  S\dX 

4tt  4lt  ' 

106  -  f  _ 

=  - A0Q  At  {  R  [  B(X)  *  T  flj  +  Rc  T  ft.dX  },  (6) 

4n  • 

by  convolving  the  instrument  function  with  the  band  brightness  function.  Here  Rs  is  the  total  band 
brightness  in  Rayleighs;  B(X)  is  normalized  band  strength;  A0  is  the  effective  area  of  the  instrument 
in  cm2;  Q  is  the  instrument  detection  quantum  efficiency;  At  is  the  integration  period  in  second  and 
Rc  is  the  continuum  background  brightness  in  Rayleigh/A.  The  instrument  detection  quantum 
efficiency  Q  is  equal  to  QkTsvs,  where  Qk  is  the  quantum  efficiency  of  the  detector  at  the  wavelength 
corresponding  to  the  k-th  channel  and  Tsv<  is  the  total  system  transmission.  If  we  consider  Qk  to  be 
independent  of  channel,  twe  only  parameters  in  equation  (6)  which  vary  w  ith  channel  number  are  Tk 
and  Qk. 


III.  Data  Inversion  and  Error  Analysis 

The  total  band  brightness  and  the  rotational  temperature  of  the  emitting  molecule  can  be 
recovered  from  the  measured  multichannel  IPD  signals.  The  instrument  function  shown  in  figure  3 
can  be  affected  by  the  thermal  state  of  the  instrument.  All  the  simulations  show  n  previously  have 
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assumed  that  the  instrument  is  at  a  reference  temperature.  As  the  temperature  changes,  the 
instrument  function  will  be  shifted,  and  the  amount  of  the  shift  AX.  is  determined  by  the  filter 
characteristics.  The  dark  counts  of  the  IPD  would  also  vary  with  the  temperature,  and  they  can 
usually  be  corrected  before  any  data  analysis  is  performed.  Here,  we  have  assumed  that  the  dark 
counts  are  very  small  and  are  excluded  in  the  present  analysis.  Therefore,  there  are  four  parameters 
which  could  be  recovered  from  the  measured  spectrum:  (1)  the  total  brightness  of  the  band  system 
R5,  (2)  the  rotational  temperature  of  the  emitter  Tr,  (3)  the  continuum  background  brightness  Rc, 
and  (4)  the  center  wavelength  shift  due  to  the  changes  of  the  thermal  state  of  the  instrument  AX. 

There  exist  many  numerical  techniques  which  can  be  used  for  our  data  analysis.  Here  we  use  a 
non-linear  least  square  fitting  technique.  In  order  to  conserve  computations,  we  linearize  the 
problem  by  expanding  equation  (6)  about  a  set  of  approximations  or  initial  guesses  (Rs(a\  Tr(a), 
Rc<a\  and  AX^)  and  keep  only  the  first  order  terms  in  the  expansion:10 

in6 

Nk  =  - —  Aq Q  At  {  R™[  B(X)* 

[B(A)*?kQkJ(a)[R[a+1)- 

R(sa)[  —  (B(X)  *?.«.)) 
dAX 

If  measured  counts  for  k-th  channel  is  expressed  as  (106AoQAt/4jt)  Nmk,  the  square  of 
deviation  from  the  theoretical  computed  signal  summing  over  all  the  anodes  can  be  calculated  by 

<5:>  =  (10*  A  ,Q  At  /4rt  r  ^  -  Nk  )2,  (7) 

k 

where  Nk  =  NJ{\Q6 A.qQA\J4k).  By  minimizing  <52>,  one  can  obtain  a  better  estimated  value  of 
Rs,  Tr,  Rc  and  AX  (denoted  by  superscript  u+1 ),  or  X(<a+1)  a  general  representation  of  the  newly 
estimated  parameters  from  a  single  matrix  multiplication  operation:13 


wol+  RrJTAdx  + 

X 

R.cx)  J  +  R(a)  f  _a_  (  B(X)  „  f^J  (a)  (T(“+1L  T(a))  + 

3T 

r 

<HI  (AXia+n  -AX(a)  )  +  (J^QjdX)  (R^+1)  -  R^a))) .  (6a) 


X"“"  =  M  1  P, 

'  u  J 


(8) 


with 
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X,  =  R 

1  s 

v  _  T(a+1) 
A2  “  ‘r 

X3  =  AX<a+1) 
X  =  R(a+1) 

4  c 


M',  -  I  A  Am  . 


pi =  i\D. 


where 


/\.  -  -!  BiX)  *  Tk^]  <“> 

\i  =  -  C' I  &<»(*>  •VVJ' 


\3  = 
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•  (B(X)*Tkflk)]W 
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-  J  TA  dX- 


*\,t  *  Rs10'  1  T,‘" I  Z*r 1  Ba> •  W  1“’ +  A*‘“’  i  <  BIX)  •  T,nk)i" 


The  degree  of  precision  obtained  depends  on  the  proximity  of  the  input  guess  values  to  the  actual 
values.  Iterations  are  performed  in  order  to  obtain  convergence  of  the  results. 


The  accuracy  of  any  measurement  is  governed  by  the  accuracy  of  the  instnr  i^nt  calibration  and 
the  counting  statistics.  Assuming  there  is  no  calibration  error,  the  standard  deviation  of  the  four 
deduced  parameters,  AX;  =  ( 6X ^ i  =  1  to  4.  can  be  simply  derived  from  the  theory  of  error 
propagation. 
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where  Nk  is  the  output  count  from  the  k-th  channel  of  the  detector  and  is  a  matrix  derived  from 
equation  (7).  We  note  that  the  random  errors  of  the  detector  can  be  described  by  a  Poisson 
distribution  and  the  standard  deviation  of  measured  signal  is  given  as  (Nk)1/2. 


IV.  An  Example:  02  (0-1)  Atmospheric  Band 

This  instrument  is  well  suited  to  measuring  the  molecular  emission  band  structure,  such  as  the 
02  Atmospheric  band,1  and  OH  Meinel  band  etc.  The  02  Atmospheric  band  is  one  of  the  strongest 
band  in  both  terrestrial  dayglow  and  nightglow  and  has  an  emission  profile  peaked  at  about  95 
km.4,14  16  It  originates  from  the  metastable  transition  between  the  excited  b'Sg+  state  and  the 
ground  state  X3Xg~  of  02  and  has  a  radiative  lifetime  of  about  12  seconds.  The  (0-0)  band  of  the 
02  Atmospheric  system  is  completely  absorbed  by  the  lower  atmosphere  and  can  only  be  observed 
from  altitude  above  the  stratosphere.  The  (0-1)  band,  however,  is  not  attenuated  thus  is  best 
suitable  for  ground-based  observation. 5'7-8  Figure  5  gives  the  emission  spectra  of  02  (0-1) 
Atmospheric  band  at  temperatures  of  150  K  and  300  K  for  an  arbitrary  band  intensity  of  1  kR.  The 
line  positions  and  strengths  are  calculated  based  upon  the  molecular  constants  of  Babcock  and 
Herzberg,17  the  intensity  measurements  of  Miller,18  and  the  Frank-Condon  factor  of  Nicholls.19 
The  band  origin  corresponds  to  a  forbidden  transition  between  two  rotationaless  states  J’  =  0  and  J" 
=  0  near  8645  A.  The  band  head  is  formed  in  the  R  branch  and  lies  on  the  short-wavelength  side  of 
the  band  origin.  The  lines  of  the  P  branch  are  located  on  the  long-wavelength  side  of  the  band 
origin.  Because  of  rotational  level  redistribution,  the  emission  spectrum  shows  different  structure 
as  temperature  changes.  Near  the  R  branch  band  head  and  the  P  branch  band  tail,  the  line  intensities 
increase  with  temperature.  This  provides  us  with  a  tool  of  monitoring  the  rotational  temperature  of 
the  excited  02(  b1£g'*' )  from  the  changes  in  band  structure  with  temperature. 

Typical  zenith  O2(0-l)  Atmospheric  band  brightness  at  night  is  approximately  500 
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Rayleighs.5,15  Figure  6  shows  the  calculated  instrument  12  channels  output  signals  at  emission 
temperatures  of  150  K  and  300  K  respectively.  The  continuum  background  emission  is  assumed  to 
be  2.0  R/A.  The  instrument  parameters  used  in  the  simulation  are  given  in  Table  1.  For  an 
integration  time  of  10  seconds,  the  maximum  signals  of  approximately  3000  counts  can  be 
obtained,  and  the  variation  in  signals  with  the  emission  temperature  is  very  clear. 

Following  equation  (8),  it  usually  takes  approximately  5  iterations  to  recover  the  total  band 
brightness,  the  temperature,  the  brightness  of  the  background  emission  and  the  filter  center 
wavelength  shift.  Table  2  illustrates  the  rate  of  convergence  for  the  spectrum  shown  in  figure  6(b), 
starting  from  a  set  of  rather  poor  guess  values. 

Figure  7(a)  presents  the  estimated  error  in  the  deduced  brightness,  ARS,  for  emission 
temperatures  of  150  K  and  300  K  as  a  function  of  brightness  Rs.  The  integration  period  used  here 
is  10  seconds,  and  the  continuum  background  is  2.0  R/A.  The  error  ARS  in  this  case  increases 
from  8  R  at  Rs  =  102  R  to  approximately  1 10  R  at  Rs  =  105  R  for  emission  temperature  of  300  K. 
This  constitutes  an  accuracy  in  the  deduced  brightness  of  less  than  \%  for  Rs  greater  than  103  R, 
and  it  improves  as  the  emission  temperature  decreases.  As  it  was  indicated  in  equation  (9),  ARS  is 
proportional  to  the  square  root  of  the  output  counts.  Therefore  ARS/(RS),/'2  approaches  a  constant 
as  Rs  increases  to  a  level  when  continuum  background  emission  is  negligible  (  Figure  7b).  Figure 
8  gives  the  error  in  the  deduced  emission  temperature.  Note  that  the  output  counts  Nk  is 
proportional  to  106AQAt/4rt,  implying  that  the  error  in  the  deduced  physical  quantities,  AX,  is 
inversely  proportional  to  the  square  root  of  the  instrument  parameters,  A0,  Q,  and  the  integration 
period  At.  By  observing  O-,  (0-1)  Atmospheric  band  system  at  night  (  ~  500  Rayleighs),  the 
instrument  described  here  is  able  to  measure  the  band  brightness  with  an  accuracy  of  ±  10 
Rayleighs  and  the  temperature  with  an  accuracy  of  ±  6  K,  for  an  integration  period  of  10  seconds, 
and  ±  4  Rayleighs  and  ±  2.6  K  for  a  integration  period  of  I  minute. 


V.  Summary 

In  this  paper,  we  introduced  a  low  resolution  spectrophotometer  with  an  interference  filter  as 
wavelength  dispersive  element  and  a  multichannel  Image  P'ane  Detector  as  light  collecting  device. 
This  instrument,  characterized  by  ns  relatively  high  throughput  and  modest  resolution,  is  suited  to 
monitoring  the  brightness  and  rotational  temperature  of  molecular  band  emission.  A  simulated 
instrument  was  presented  here  as  an  example  to  demonstrate  its  capability  in  measuring  the  O,  (0-1) 
Atmospheric  band  emission  in  the  mghtglow.  For  an  integration  period  of  10  seconds,  this 
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instrument  is  able  to  measure  temperature  with  an  accuracy  of  ±  6  K  and  brightness  with  an 
accuracy  of  ±  10  Rayleighs  for  a  typical  zenith  emission  brightness  of  500  Rayleighs. 
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Table  1.  Parameters  of  Simulated  Instrument 


Interference  Filter 

center  wavelength,  Xo 

full  width  at  half  height  (FWHH) 

maximum  transmittance 

Objective  Lens 

focal  length 


1PD 

anode  structure 

area  of  anode 

radius  of  anode  ring 

quantum  efficiency  at  8730  ±  50 A,  8k 

Aperture 

effective  area,  Ao 


8730  A 
-12  A 

0.5 


2.27  cm 


12  concentric  ring, 
equal  area  anodes 
0.076  cm2 

(see  Killeen  era/.  1983) 
0.005 


1  cm2 


Total  System  Transmission,  Tsys 


0.4 


I  able  2.  An  Example  of  Data  Analysis  Output 
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FIGURE  CAPTIONS 


Fig  1:  Schematics  of  the  Image  Plane  Detector  Spectrophotometer. 

Fig  2:  The  transmission  function  of  a  simulated  double  cavity  interference  filter  with  peak 
transmission  normalized  to  1 .0. 

Fig  3:  The  filter  transmission  functions  for  12  light  incident  angles:  8k  =  (k-l)1/29c,  k  =  1,  2, .... 
12  and  9C  =  4.4°. 

Fig  4:  Instrument  transmission  function  Tk,  multiplying  the  field  of  view  (k  =  1,  2 . 12). 

Fig  5:  O2(0-l)  Atmospheric  band  emission  spectra  with  total  band  brightness  of  lkR  at  a 

temperature  of  (a)  150  K  and  (b)  300  K. 

Fig  6:  The  calculated  12  channel  LPD  Spectrophotometer  signals  for  a  (0-1)  band  source 
brightness  of  500  Rayleighs  at  an  emission  temperature  of  (a)  150  K  and  (b)  300  K 
(integration  time  =  10  seconds). 

Fig  7:  (a)  Absolute  error  ARS,  and  (b)  normalized  error  ARS  /(  R$)1/2  in  the  deduced  brightness 

for  an  integration  period  of  10  seconds. 

Fig  8:  Absolute  error  in  the  deduced  temperature  for  an  integration  period  of  10  seconds. 
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